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ABSTRACT 


A research program is being conducted to study the mechanisms controlling the spread 
of fire in turbulent forced gas flows moving in the direction of flame propagation. The 
research tasks completed during this reporting period include experimental studies of the 
effect of flow turbulence on the rate of concurrent flame spread, and of mass burning when 
the combustible material is in a floor configuration. The results of the experiments with thick 
PMMA sheets show that flow turbulence affects significantly the flame spread and mass burn- 
ing processes. As the turbulence intensity is increased, the flame spread rate decreases because 
the flame length sharply decreases, and the mass burning rate increases because the surface 
heat flux increases. Currently underway is a complementary study with the combustible 
material placed in a ceiling geometry to observe the effect of buoyancy on the flame spread 
and mass burning processes. An additional task also completed during this period is a review 
of the processes of ignition and flame spread of solid combustibles. The review presents a 
novel approach to the analysis and prediction of both processes based in their close inter- 


relationship. 
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I. INTRODUCTION 


During this reporting period we have been conducting experiments aimed to study the 
effect of flow turbulence on the spread of fire over the surface of combustible materials. The 
study is of practical interest because in most realistic situations the spread of fire occurs under 
turbulent flow conditions, and there is a need to understand and determine how turbulence 
affects the mechanisms that control the spread process. Furthermore, the experiments are 
yielding results that are potentially important not only in flame spread modeling and predic- 
tion, but in other aspects of fire development and testing such as burning rates, flame length 
and soot emission. Particularly interesting are the results obtained for the concurrent mode of 
flame spread, where significant observations have been made concerning the influence of flow 
turbulence intensity on the rate of flame spread, the mass burning rate, and the flame length. 
The results are specially important because the concurrent mode of flame spread, being the 
fastest and most hazardous among the different modes of flame spread, has a great importance 


in the fire safety field. 


A series of experiments to measure the rate of mass burning and flame spread over thick 
PMMA samples as a function of the velocity and turbulent intensity of the air flow have been 
completed during the present grant period. The solid fuel sample is in a floor type 
configuration, and the flow turbulence is generated in the free stream. Currently we are con- 
ducting similar experiments but with low velocity flows and with the PMMA samples placed 
both on a floor and a ceiling configurations to observe the effect of buoyancy on both tur- 
bulent mass and flame spread rates. A summary of the research progress made during the 
current grant period is presented below, and a more detailed presentation is given in the 


appendices. 


Il. RESEARCH PROGRESS 
II.1 Concurrent Turbulent Flame Spread 


An experimental study has been conducted of the influence on the spread of flames over 
the surface of PMMA sheets of the turbulence intensity of an air flow in the same direction of 
flame propagation. The experiments are carried out in a facility specifically designed to per- 
form flame spread experiments consisting of a small scale combustion wind tunnel and sup- 
porting instrumentation based primanily on optical and thermocouple measuring methods. The 
turbulence intensity in the tunnel test section is varied by means of grids and perforated plates 
of different sizes placed at the exit of the tunnel converging nozzle. The velocity and tur- 
bulence intensity of the flow field are measured with a one component Laser Doppler Velo- 
cimeter. The flame spread rate is measured from the surface temperature history as given by 
thermocouples embedded at fixed intervals along the fuel surface. A Schlieren system is used 
to provide qualitative information about the effect of the flow turbulence on the flame and 
thermal layer structure. The fuel specimens are 6 in. wide, 12 in. long and 0.5 in. thick 
PMMA sheets and the oxidizing gas is air. Tests are conducted for flow velocity ranging 


from | to 4 m/sec with turbulence intensity varying from 1 to 20%. 


The measurements of the flame spread rate over PMMA sheets as a function of the con- 
current air flow velocity are shown in Fig. | for several values of the turbulence intensity. It 
is seen that the spread rate increases approximately linearly with the flow velocity, that it 
decreases as the turbulence intensity increases, and that this last effect is more pronounced for 
larger flow velocities. These results are very interesing and somewhat surprising since this 
mode of flame spread is controlled by heat transfer from the flame to the fuel, and it is well 
known that turbulent boundary layer heat transfer is larger than the laminar flow one. The 


results, which appear to be due to a strong effect of the turbulence intensity on the flame 
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length, are very important not only because they introduce new aspects about the flame spread 
process not previously predicted, but because it may have significant influence in the applica- 


tion of flame spread formulas in models of room fire development. 


The mechanisms by which turbulence affects the flame spread rate can be inferred from 
the theoretical analysis of the spread process. A simplified heat transfer model of the flame 
spread provides the following expression for the rate of spread (Saito, L., Quintiere, J. and 


Williams, F.A., 1st Int. Symp. Fire Safety Sci. 75, 1986). 


ve = q? Ip(kpe (T, - T;))? (1) 


where q is the surface heat flux, |; the flame length, kpc are the thermal properties of the solid 
and T,, and T; the solid pyrolysis and initial temperatures respectively. The flow velocity and 
turbulence intensity can affect both q and |, and through them the flame spread rate. Thus, it 
is important to determine how turbulence affects these parameters. In the work performed to 
date, we have not measured these effects directly. However, it is possible to deduce them 
approximately from the surface temperature histories. The flame length is determined with the 
spread rate and the time required for the surface temperature at a thermocouple position to 
rise from ambient to the pyrolysis value. The surface heat flux is calculated from the time 


variation of the surface temperature by assuming that the fuel behaves as a semi-infinite solid. 


The calculated variation of the flame length with the flow velocity and turbulence inten- 
sity 1S presented in Fig. 2. It is seen that the dependence of the flame length on the flow 
velocity is different depending on the turbulence intensity. For low turbulence intensities, the 
flame length increases with the flow velocity and for large turbulence intensities, it decreases. 
The calculated variation of the total heat flux with the turbulence intensity is shown in Fig. 3. 
It is seen that the heat flux is only weakly affected by the flow turbulence, decreasing slightly 


with the turbulence intensity for large flow velocities and increasing slightly for low flow 
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velocities. All of the above results have been combined in Fig. 4 where the non-dimensional 
flame spread rate deduced from Eq. (1) is plotted versus the flow turbulence intensity. It is 
seen that the non-dimensinal spread rate is independent of the flow velocity and turbulence 
intensity, which shows the validity of Eq. (1) and verifies that the effect of the flow tur- 
bulence on the flame spread rate takes place primarily through the flame length and heat flux. 
This again is a significant result since the prediction of flame length as a function of the prob- 


lem parameters is an important factor in the development of room fires. 


A detailed description of this work is given in Appendix A. The Appendix is a paper 


that has been submitted for presentation at the 23rd Int. Symp. on Combustion. 


11.2 Turbulent Surface Burning 


Measurements have been made of the PMMA burning rate dependence on the velocity 
and turbulence intensity of the air flow. The objective of the measurements is to provide 
additional information about the mechanisms that control the observed effect of flow tur- 
bulence on the flame spread rate. As it was reported above, the concurrent flame spread rate 
decreases as the turbulence intensity is increased due primarily to the shortening of the flame 
length with the turbulence intensity. Since the flame length is the result of the combustion of 
the excess pyrolyzate produced at the pyrolysis region, it is important to determine how the 


regression rate, and consequently the pyrolysis rate, depends on the flow turbulence. 


The measured regression rate of the PMMA sheet for a flow velocity of 2 m/sec is 
shown in Fig. 5 as a function of the distance from the specimen upstream leading edge, with 
the turbulence intensity as parameter. The experimental data for other flow velocities show 
very Similar trends and are not presented here. Each data point is an average of two tests 


under the same flow conditons. From Fig. 5, it is seen that. the regression rate becomes 
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this parameter is plotted versus Re?*(u’/U), a correlation between these two non-dimensional 
parameters is observed. Though this correlation can be qualitatively explained in terms of the 
turbulent heat transfer mechanisms that control the burning process, the reason why the flow 


problem parameter take the form indicated in Fig. 6 is yet to be explained. 


A detailed description of this work is given in Appendix B, which presents a paper 
prepared for presentation at the Third International Symposium on Fire Safety Science, Edin- 


burgh, Scotland, July 1991. 


II.3 The Solid Phase: On Ignition and Flame Spread 


This work is basically a review of ignition and flame spread in solid combustibles, which 
also includes a novel analysis that provides general engineering formulas for ignition delay 
times and flame spread rates. The work has been prepared for its publication in a Combustion 


Treatise on Fire edited by Mr. Geoff Cox, and that will be published by Pergamon Press. 


The novel contribution of the work is the realization that the flame spread rate can be 
expressed as the ratio of the length of the solid heated region ahead of the pyrolysis front to 
the solid ignition delay time. Thus, if an analysis of ignition is developed for a given 
geometry and flow configuration, the corresponding flame spread analysis can be derived from 
the ignition one by simply calculating the solid heated length. In the work, an analysis of 
solid fuel ignition in a flat plate flow is developed that incorporates the solid heating and gas 
phase reaction controlling mechanisms, and that provides a general formula for the ignition 
delay time. This analysis is then extended to obtain formulas for the flame spread rate that 


describe well the experimental observations. This work is presented in Appendix C. 
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ABSTRACT 


The rate of flame spread over the surface of thick PMMA sheets has been measured as a 
funcuon of the velocity and turbulence intensity of a forced air flow concurrent with the direc- 
tion of flame propagation. It is shown that the flow turbulence has a strong influence on the 
flame spread process. For all turbulence intensities, the flame spread rate increases monotoni- 
cally with the flow velocity, although the rate of increase is smaller for larger turbulence intensi- 
ties. For a given flow velocity, the spread rate decreases as the turbulence intensity is increased, 
the effect being more pronounced for larger flow velocities. These effects appear to be due to a 
strong influence of the turbulence intensity on the flame length, which decreases as the tur- 
bulence intensity is increased, thus reducing the net heat flux from the flame to the unburnt com- 
busuble. The flame spread rate measurements are correlated well with an expression for the 
pyrolysis front rate of spread that is deduced from a simple solid phase energy analysis, which 
verifies the above observations and confirm the conclusions from previous works that heat 
transfer from the flame to the solid fuel is the controlling mechanism in this mode of flame 
spread. The results are significant because the prediction of the flame length and spread rate as a 
function of the problem parameters are important factors in the development of room fire 
models, and the establishment of material flammability tests. 


INTRODUCTION 


In the concurrent mode of flame spread, the oxidizer gas flows in the spread direction thus 
assisting the spread of the flame. The oxidizer flow can be induced naturally by the flame 
(upward flame spread) or, as considered in this case, induced by externally applied pressure gra- 
dients (forced flow flame spread). Among the different modes of flame spread, this is the fastest 
and most hazardous because the gas pushes the flame ahead of the burning region, which 
enhances the heat transfer from the flame to the unburnt material and consequently the spread of 
the flame. This type of flame spread is often encountered during the spread of fire, where the 
flow conditions are generally turbulent. Thus, the study of concurrent turbulent flame spread is 


of interest in the fire safety field. 


The concurrent, or flow assisted, flame spread is characterized by the following mechan- 
isms (see insert in Fig. 1). Upstream from the pyrolysis front heat transfer from the flame to the 
solid causes it to pyrolyze and establishes an approximately constant temperature through a bal- 
ance between heat transfer and fuel vaporization. The vaporized fuel is diffused and convected 
outward and forward reacting with the ambient oxidizer and establishing a diffusion flame that 
lies in the boundary layer next to the solid surface. The fuel vapor that is not consumed in the 
upstream flame is convected downstream from the pyrolysis front where it keeps reacung with 
the oxidizer, extending the diffusion flame downstream. Heat transfer from the flame to the 
unburnt fuel downstream from the pyrolysis front causes the solid temperature to increase until it 
reaches its pyrolysis temperature at which point the solid starts to pyrolyze. The onset of fuel 
pyrolysis determines the progress of the pyrolysis front and consequently the flame spread rate. 
Thus, the flame spread process is controlled by the transfer of heat from the flame to the solid 
fuel, and by the physico-chemical mechanisms that generate the diffusion flame. Since flow tur- 
bulence affects both controlling mechanisms, it is important to determine how they are affected 


by turbulence in order to predict the spread of flames under realistic conditions. 


The only works published to date on turbulent flame spread are those of Orloff et al.! for 


upward turbulent flame spread over large PMMA sheets, and of Saito and co-workers? for 


upward turbulent spread of flames over PMMA and wood sheets. In these works, the gas flow is 
buoyancy induced, and the turbulence is naturally generated depending on the sample scale. The 
studies concentrate on measuring the rate of flame spread dependence on scale! or the magnitude 
of an imposed radiant flux*. With the exception of the study of Zhou et al.4 for opposed flow 
flame spread, no systematic studies have been performed to date to observe the effect of flow 
turbulence intensity on the rate at which flames spread over the surface of a solid fuel. In the 
present work, a parametric study is carned out of the spread of flames in a concurrent turbulent 
forced flow with the objective of determining what is the effect of turbulence intensity on the 
rate of flame spread, and to infer through these measurements how turbulence affects the dif- 


ferent mechanisms controlling the spread of the flame. 


EXPERIMENT 


A schematic diagram of the expenmental facility is shown in Fig. 1. It consists of a labora- 
tory scale combustion tunnel, designed to conduct combustion experiments under vaned flow 
conditions, and the supporting instrumentation, based primamily on optical and thermocouple 
measuring methods. The tunnel has a 0.89 m long settling chamber with a rectangular cross sec- 
tion 0.3 m by 0.18 m, which supplies a gas flow to the tunnel test section through a converging 
nozzle with a reduction area of 5.5 to 1. The test section 1s 0.6 m long and has a rectangular 
cross section 130 mm wide and 76 mm high. The side walls are made of 6 mm Pyrex glass to 
enable visual observation and optical diagnostic access, and the floor and ceiling sides are made 
of 25.4 mm thick Marinite sheet. The tunnel is mounted horizontally on a three axis motonzed 


positioning table. 


The gas flow in the tunnel is supplied from centralized compressed air and the flow rates 
controlled and measured with critical nozzles. Turbulence is generated by means of grids, or 
perforated plates, placed at the exit of the tunnel converging nozzle. A prescribed turbulence 
intensity is obtained through a combination of flow velocity and plate blockage ratio. This 


means of turbulence generation. however, cannot provide high turbulence intensity at very low 


flow velocities. It also causes a change in turbulence scale as well as intensity, which will affect 
the turbulent mixing. Information about the distribution of turbulence intensity through the test 
section is given in Ref. [4] and will not be repeated here. It is characterized by an initial region, 
approximately 50 mm where the turbulence decays sharply, followed by a 0.2 m region of rela- 
tively constant turbulence (where the measurements are primarily made), and a final region 
where the turbulence decays again due to entrainment effects at the test section exit. The tur- 
bulence intensity distribution normal to the fuel surface 1s similar to that of the velocity, constant 
in the center region and decreasing toward fuel surface. In this work, all the data are correlated 
in terms of the constant free stream velocity and turbulence intensity. The gas velocity and tur- 
bulence intensity are measured with a one component Laser Doppler Velocimeter (LDV) operat- 
ing in the forward scattered mode. The laser source is a Spectra Physics Ar-Ion laser with a 2 W 
total ourput power. The signal is processed with a TSI frequency counter and transferred to an 
IBM/AT micro-computer. The experimental installanon also includes a Schlieren system with a 
0.5 m diameter collimated light beam that is used to provide qualitative information about the 
effect of the flow turbulence on the flame and thermal layer structure. The system consists of a 
Tungsten lamp light source, beam expanding and converging lenses and two 0.5 m parabolic 
murrors with a4 m focal length. The Schlieren images are recorded with still photography and 
high resolution video. In addition the installation contains a network of thermocouples for gas 


and solid phase temperature measurements. 


The fuel specimens tested in this work, are 12.7 cm thick PMMA (Rohm and Hass, Plexig- 
las G) sheets, 76 mm wide by 0.3 m long, which are considered as representative of a thermally 
thick solid. They are mounted flush in the tunnel Marinite floor with their upstream edge placed 
15 mm from the tunnel convergent nozzle exit. The flame spread process is initiated by igniting 
the top upstream edge of the sheet with an electrically heated nichrome wire placed in close con- 
tact with the PMMA. The flame spread rate is measured from the surface temperature histories 
as given by thermocouple placed at fixed intervals along the fuel surface. Eight Ch.-Al. thermo- 
couples 0.127 mm diameter are embedded on the PMMA with their beads flush with the PMMA 


Surface at distance 32 mm apart. The thermocouples output is processed in the micro-computer. 
With the surface temperature histories, the rate of spread of the pyrolysis front is calculated from 
the time lapse of pyrolysis arrival to two consecutive thermocouples and the known distance 
between the thermocouples. As it is explained later, they are also used to calculate flame and 


pyrolysis lengths and surface heat fluxes. 


RESULTS AND DISCUSSION 


In order to better understand the characteristics of the experimental results, it is convenient 
to bnefly examine the mechanisms determining the spread of the flame. Previous experimental 
and theorencal work on the concurrent mode of flame spread indicate that heat transfer from the 
flame to the solid fuel is the dominant controlling mechanism>®. In fact, a simple energy 
analysis applied to a control volume in the unburnt solid downstream from the pyrolysis front 
provides an expression for the flame spread rate’ that seems to describe (by convection and radi- 
ation) the spread process well. Assuming that the solid behaves as semi-infinite, that the heat 
flux from the flame to the solid, q¢ is constant throughout the downstream flame length, /;, and 
that the solid pyrolyzes (and ignites) when its surface temperature reaches the pyrolysis tempera- 
ture, Tp, the following expression is obtained for the spread rate2’, Wis 


Viz 4( gel d/2)2(1e/lp) 


PR kpe (Tp —T)) ms 


where /, is the pyrolysis length, kpc are the thermal properties of the solid (thermal inertia) and 
T; the initial solid temperature. From Eq. (1) it is seen that the flow turbulence can affect the 
flame spread rate through the heat flux and the flame length primarily, although other mechan- 


isms are also possible. 


The measurements of the flame spread rate over PMMA sheets as a function of the con- 
current air flow free stream velocity are shown in Fig. 2 for several values of the free stream tur- 
bulence intensity. The spread rate is an average of the values deduced from consecutive thermo- 


couples throughout the specimen length and from three different tests. The standard deviation 1s, 
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in most cases, of the order of 7%. No increasing or decreasing trend could be eee from the 
local measurement and consequently, the spread rate was considered as constant throughout the 
specimen. Measurements for flow velocities smaller than 1 m/sec were not made because buoy- 
ancy had too strong an influence in the flame spread process. At velocities larger than 4 m/sec, 
signs of exuncton appeared at the upstream leading edge. It is seen that the spread rate 
increases monotonically with the flow velociry, but that the rate of increase is smaller as the tur- 
bulence intensity increases. The monotonic increase of the spread rate with the flow velocity is 
predicted by most of the present models of concurrent flame spread®, and, as it is explained later. 
is consistent with heat transfer processes in forced flow flat plate boundary layer. As the flow 
velocity increases, the thickness of the boundary layer, and consequently the flame stand-off dis- 
tance, decreases which results in a larger surface heat flux and in turn, on a larger flame spread 
rate. The dependence of the concurrent flame spread rate on the turbulence intensity is presented 
in Fig. 3. It is seen that the spread rate decreases as the turbulence intensity increases, and that 
the effect is more pronounced for larger flow velocities. These results are very interesting and 
somewhat surprising since this mode of flame spread is controlled by heat transfer from the 
flame to the fuel, and it is well known that turbulent boundary layer heat transfer is larger than 
the laminar flow one. The results are important not only because they introduce new aspects of 
the flame spread process not previously predicted, but because they may have significant 


influence in the application of flame spread formulas in models of room fire development. 


In order to understand the above results, it is important to determine how the flow velocity 
and turbulence affect the flame spread rate through both q¢ and /;. In the work performed here, 
we have not measured these effects directly, but we have deduced them approximately from the 
surface temperature histories. The flame length can be determined from the tme required for the 
surface temperature at a thermocouple position to nse from ambient to the pyrolysis value in 
conjunction with the flame spread rate. It should be noted that the resulting length is not neces- 
sarily the actual flame length, but the length of the elevated temperature region ahead of the 


pyrolysis front. For simplicity this length is called here, flame length. Using this method, the 
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vanianion of the flame length to pyrolysis length ratio with the flow velocity and turbulence inten- 
Sity were determined and are presented in Figs. 4 and 5 respectively. The flame and pyrolysis 
length values are calculated from data provided by eight thermocouples placed at equal distances 
along the PMMA surface. The ratio of both lengths was found to be approximately constant 
along the PMMA surface, and the value presented in figures correspond to an average of all the 
calculated ranos. From the results of Fig. 4, it is seen that the flame length to pyrolysis length 
ratlo increases approximately linearly with the flow velocity although the rate of increase is 
smaller for high turbulence intensity. From the results of Fig. 5, it is seen that for all flow velo- 


cities, the flame length to pyrolysis length rano decreases as the turbulence intensity increases. 


Schlieren images taken during the spread of the flame were also used to observe the effect 
of flow turbulence on the flame length and the thermal boundary layer. The images show that 
the thermal boundary layer has the characteristics of a flat plate forced flow one with its thick- 
ness increasing with the distance from the upstream leading edge. They also show the strong 
inducement of cold air into the flame zone by the flow turbulence, which results in the eventual 
shortening of the flame. The flame shortening can be the result of either quenching or enhance- 
ment of the combustion reaction. In the former case, the strong inducement of cold air into the 
reaction zone, and the interaction of the cold solid wall with the flame, could quench the reaction 
downstream where the fuel is less abundant and consequently, shorten the flame. In the latter 
case, the increased mixing of the fuel vapor and oxidizer could enhance the reaction consuming 
the fuel sooner, which also would shorten the flame. At this point, we do not have enough infor- 


mation to determine which one of the two mechanisms is the dominant one. 


Another important parameter in the problem is the surface heat flux in the unburnt solid 
downstream from the pyrolysis front. In this work, we have calculated it from the variation with 
time of the surface temperature by assuming that the fuel behaves as a semi-infinite solid 
exposed to a constant surface heat flux. This method is applied to the data provided by each one 
of the surface thermocouples to determine the variations of the heat flux along the surface and 


obtain the corresponding correlation. It was found that the product qe /;!/? was approximately 
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constant along the surface, in qualitative agreement with forced flow boundary layer heat 
transfer predictions. The calculated variation of the heat flux with the flow velocity is shown in 
Fig. 6 for several turbulent intensities. The heat flux data is an average value of those obtained 
through the fuel surface and corresponds approximately to that received by a surface element 
that is midway between the flame tip and the pyrolysis front. The results show that the surface 
heat flux increases monotonically with the flow velocity although the increase is less marked for 
higher turbulent intensities. For a fixed flow velocity the surface heat flux increases with the tur- 


bulence intensity. 


In order to verify the predicted capabilities of Eq. (1), all of the results from Figs. 2, 4 and 6 
have been combined to produce Fig. 7, where the non-dimensional flame spread deduced from 
Eq. (1) is plotted versus the flow turbulence intensity for all the flow velocities tested. It is seen 
that the correlation of the measurements 1s remarkably good particularly considering the 
sumplified character of Eq. (1) and the diversity of the measurements. The deviation from unity 
in the non-dimensional flame spread rate and the scatter on the correlation is attnbutable to the 
selection of the T, and kpc PMMA (taken as constant with T, = 390°C, k = 1.99x10-? J/smK, 

p = 1.99 kg/m}, c = 1.46 kJ/kg K) and the approximate determination of the flame lengths and 
surface heat fluxes. A value of T, = 300°C would result in a unity non-dimensinal flame spread. 
We feel, however, that such value is too low to represent the PMMA pyrolysis temperature. The 
results of Fig. 7 in conjunction with those of Figs. 4 and 6 corroborate the above arguments that 
the shortening of the flame length with the turbulence intensity is the primary reason for the 


observed reduction of the flame spread rate as the turbulence intensity increases. 


CONCLUDING REMARKS 


The results of this study show that flow turbulence has a strong influence on the flame 
length, surface heat flux, and on the concurrent flame spread rate. Both the observed shortening 
of the flame length and decrease of the flame spread rate as the turbulent intensity is increased 


are significant new observations that may have an important impact in fire modeling. since the 


22 


prediction of flame lengths and rates of flame spread as a function of the environmental condi- 
tions are important factors in the development of room fire models and the establishment of 


material flammability tests. 


The limitations of this experiment; specimen size, measurement difficulties, indirect data, 
etc., did not allow us to extract more accurate information about the problem. For example, a 


convection based heat transfer analysis applied to Eq. (1) predicts 


Ve (Tr- Tp)? le O 
(ie Dosueny dp ies c 


where the surface heat flux has been taken as qe =c Re® Pr!4 (T; - T)/lp according to boundary 
layer theory. In the above expressions, T; is the flame temperature, U is the free stream flow 
velocity, Re and Pr are the Reynolds and Prandtl numbers, a is an exponent that may vary for 
laminar and turbulent flows and C is a constant. Thus, if ///, remains constant along the fuel 
Surface, Eq. (2) predicts a constant spread rate that increases linearly with U.. for laminar flow, 
and one that accelerates as it spreads over the solid surface and has a power law dependence on 
U.. for turbulent flow. Our present data do not show an accelerating trend in turbulent flow 
perhaps because the data is not detailed enough due to the above indicated experimental limita- 
tions. However, it should be noted that in this work, the free stream velocity 1s relatively small 
and that turbulence is generated by introducing vorticity in the downstream free stream flow. 
Thus, the classical heat transfer approach of correlating Nusselt and Reynolds numbers may not 
totally apply here unless turbulent intensity related scale is introduced in the problem. Another 
aspect of this study that needs further work is the determination of how turbulence affects the 
flame length. In this regard, the measurement of the combustion efficiency by exhaust gas 
analysis could provide information about whether enhanced combustion or quenching cause the 
shortening of the flame with the turbulence intensity. Given the potential impact of the present 
results on the theoretical prediction of the spread of flames, it is important to continue the work 


to study these and other aspects of the problem. 
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LIST OF FIGURES 


Schemauc of experimental facility. Inset: schematic of concurrent flame spread. 


Dependence of the flame spread rate on the free stream flow velocity for several values of 
the turbulence intensity. 


Variation of the flame spread rate with the turbulence intensity for several values of the free 
stream flow velocity. 


Dependence of the downstream flame length to pyrolysis length ratio on the free steam 
flow velocity for several values of the turbulence intensity. 


Variation of the downstream flame length to pyrolysis length ratio with the turbulence 
intensity for several values of the free stream flow velocity. 


Vanation of the surface heat flux at a position midway along the flame with the free stream 
flow velocity for several values of the turbulence intensity. 


Correlation of the flame spread data in terms of a non-dimensional flame spread rate 
deduced from Eq. (1). 
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ABSTRACT 


Experiments have been conducted to study the effect of flow velocity and grid-generated 
turbulence intensity on the surface regression rate of thick PMMA sheets burning in a forced air 
flow. All the tests are carried out in a laboratory-scale combustion tunnel with flow velocities 
ranging from | m/s to 4 m/s and turbulence intensities from 1% to 20%. It is found that for all 
turbulence intensities, the regression rate decreases with the distance from the fuel sheet 
upstream edge, and increases with the flow velocity, in agreement with boundary layer analyses 
of the process. It is also found that flow turbulence has a strong influence on the surface regres- 
sion rate due primarily to an enhancement of the heat transfer from the flame to the fuel that is 
caused by turbulence generated flame fluctuations. It is shown that the experimental data can be 
correlated in terms of a non-dimensional mass burning rate, m” Lx/(A (T;— Tp)) that is approxi- 
mately linearly proportional to the turbulent flow parameter ((u’/u).. Re?-8)95. The good correla- 
tion of the data in terms of these parameters indicate the potential predictive capabilities of 
theoretical models of the process that incorporate the flow turbulence through an eddy viscosity, 
and thermal and mass diffusivities. Explicit expressions for the laminar and turbulent mass 


burning rates in terms of the solid fuel and flow properties are reported. 
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NOMENCLATURE 


a Constant 

Cp Gas specific heat 

f Normalized stream function 

Ie Heat of pyrolysis 

M Molecular weight 

me Mass burning rate 

Q Heat of combustion 

qe Convective surface heat flux 

ape Radiative surface heat flux 

Rey Reynolds number, u.. x/v 

r Surface regression rate 

ae Temperature 

u Velocity parallel to the fuel surface 
(u’/U).. Turbulence intensity at the free stream 
Vv Velocity normal to the fuel surface 

x Coordinate px. -allel to the fuel surface 
VE Mass fraction 

y Coordinate normal to the fuel surface 
Z Mixture fraction, (B — B..)/(Bp — Bu.) 


Greek Symbols 


a Thermal diffusivity 

Br Species Shvab-Zeldovich variable, Ye/Vp Mg — Yo/Vo My 

Br Temperature-species Shvab-Zeldovich variable, c,T/Q — Yo/Vo Mo 
€ Eddy kinematic viscosity and diffusivity, a(u’/u).. u. y 

n Similarity variable, Vu_/vx 


Thermal conductivity 
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Vv Kinematic viscosity, also stoichiometric coefficient 


Ne) Density 

y Stream function, Vv un x f(1) 
Subscmpts 

S Fuel 

f Flame 

) Oxygen 

p Pyrolysis 

S Solid 

oo Free stream 

- Gas phase 
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INTRODUCTION 


The surface regression, or mass burning rate of a solid combustible burning in a convective 
oxidizing gas flow is an important subject in problems related to the development of fires, 
material processing and rocket propulsion. Specifically in a fire, it determines such aspects of 
the problem as flame heights, rate of heat release, and smoke and CO emission, among others, 
which in turn determine the growth rate of the fire and its potential hazard. A configuration of 
this combustion process encountered in practical situations and amiable to theoretical modeling, 
is that of a flat solid fuel surface burning in an oxidizer gas flowing parallel to the solid surface. 
The resulting combustion problem is that of a diffusion flame established in the boundary layer 
adjacent to the fuel surface. The solid pyrolysis 1s caused and sustained by heat transfer from the 
diffusion flame to the solid surface. The characteristics of the diffusion flame and of the heat 
transferred from the flame to the surface, and consequently of the mass burning rate, depend 
strongly on the characteristics of the gaseous flow in addition to the solid properties. Since most 
practical combustion processes occur under turbulent flow conditions, it is of interest to investi- 
gate how the flow turbulence affects the heat transfer and chemical kinetic mechanisms that con- 


trol the mass burning process. 


A number of theoretical and experimental studies have been conducted of the laminar and 
turbulent burning of a flat fuel surface with a focus on fire applications. Emmons! in his analysis 
of the combustion of a flat liquid surface in a forced, laminar oxidizing flow established the 
foundations for the theoretical modeling of this problem. Using the boundary layer and flame 
sheet approximations he obtained explicit formulas for the mass burning rate in terms of the gas 
flow parameters and fuel properties. Kosdom et. al.* and Kim et. al.>, applied the same metho- 
dology to develop analyses of the laminar, free convective burning of a vertical fuel surface, and 
conducted experiments that verified the mass burning rate predictions of the model. Mao et. al.4 
have developed a mixed convection formulation of the problem that gives the mass burning rate 
for the whole range of flow conditions from forced to free and for surfaces with different onenta- 


tions. Sibulkin et. al.5 have incorporated finite rate chemistry to the model and Mao et. al.® have 


38 


developed a numenical analysis of the problem that removes both the boundary layer and flame 
sheet approximations. Further details on the problem can be found in the reviews of Pagni’ and 
Sibulkin’. The theoretical formulation of the turbulent free convective burning of a vertical fuel 
surface has been developed by Kennedy and Plumb? and Tamanini!® using infinite chemistry 
and and Sibulkin® using finite rate kinetics. Experimental measurement and analysis can also be 
found from works done by Ahmad and Faeth!!, De Ris and Orloff!? and Orloff et. al.!3.14 
among many others. These last studies focused on the mass burning rate of solid combustibles 
in a buoyantly-generated turbulent flow. The non-dimensional burning rate was generally corre- 
lated with the Rayleigh Number, and in general, the experimental data agreed well with the 


theoretical predictions. 


Another area related to the present work is the study of hybrid rocket fuel combustion. 
Notable among the different works published to date on the turbulent flow mass burning problem 
are those of Marxman and Gilbert!5, Marxman!6 and Paul et. al.!7 among others. The results of 
these works are generally presented in the form of relations between the regression rate and the 
mass transfer number B. There are some similarities between these works and the present work, 
although these works are mainly concermed with hybrid rocket propellant burning in high speed 
flow and the present work is concerned primarily with the fire aspects of the process where low 


velocity flows prevail. 


The turbulence in the studies mentioned above is either shear-generated or buoyantly- 
generated, and the flow structure is usually very complex because turbulence generation and dis- 
sipation occur simultaneously. The other alternative to study flow turbulent effects is to intro- 
duce the turbulence artificially into the flow by means of a grid or perforated plate placed per- 
pendicular to the flow direction. In grid-generated turbulent flow, the mean velocity gradient is 
very small so that the flow structure can be considered as isotropic and homogeneous, and, there- 
fore, also relatively easier to model theoretically. Although considerable experimental and 
numerical works have been done on both premixed and diffusion free flames in shear-generated 


turbulent flows!8, the turbulent flow structure is so complex in shear-generated turbulent flow 
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that the interaction between the flow and the flame is still far from clear. Therefore, it is useful 
to use simple turbulent flow such as grid flow to determine the relations between the flow and the 
flame. Even so, grid-generated turbulent diffusion flames have rarely been investigated, as indi- 


cated by Manako et. al.!9. 


In the present work a systematic study has been carried out of the burning of a solid fuel in 
a forced, grid-generated turbulent air flow with the objective of experimentally determining the 
effect of the external flow velocity and turbulence intensity on the mass burning rate, and of 


finding proper correlations between the burning rate and the flow parameters. 


EXPERIMENT 


The experimental apparatus is shown schematically in Fig. 1. It consists of a laboratory 
scale wind tunnel designed to conduct combustion experiments under vaned flow conditions, 
and the supporting instrumentation. The wind tunnel has a 0.89 m long settling chamber with a 
rectangular cross section 0.31 m by 0.18 m, which supplies gas flow to the tunnel test section 
through a converging nozzle with a reduction area of 5.5 to 1. The test section is 0.61 m long 
and has rectangular cross section 0.122 m wide and 76 mm high. The side walls of the test sec- 
tion are made of 6 mm thick Pyrex glass for visual observation and optical diagnostic access, 
and the floor and ceiling are made of 55 mm thick Marinite sheet. The combustion tunnel is 


positioned horizontally on a three axis positioning table. 


The air flow in the test section is supplied from a centralized gas compressor and the flow 
rates controlled with critical flow nozzles. Turbulence is generated by means of grids, or per- 
forated placed at the exit of the tunnel converging nozzle. A prescribed turbulence intensity is 
obtained through a combination of flow velocity and plate blockage ratio. A detailed descnption 
of the velocity and turbulence profiles through the test section is given in Ref. [20]. The air flow 
velocity and turbulence intensity are measured with a one-component Laser Doppler Velocime- 
ter operating in the forward scattered mode. The experimental installation also includes a 
Schlieren system with a 0.45 m diameter collimated light beam and an array of eight thermocou- 


ples that is used to measure and monitor the solid combustible surface temperature. 
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The fuel specimens used in this work are 12.7 mm thick PMMA (polymethyl! methacrylate) 
(Roam and Hass, Plexiglas G) sheets, 76 mm wide by 0.37 m long. PMMA was selected 
because of its non-charring burning. The fuel samples are mounted flush in the tunnel Marinite 
floor with their upstream edge placed 0.1 m from the tunnel convergent nozzle exit. The PMMA 
burning is initiated by igniting the upstream edge of the sheet with an electrically heated 
nichrome wire. The flame is left to spread over the solid surface and the arrival of the pyrolysis 
front at a given thermocouple is recorded to mark the initiation of the mass burning process at 
that particular location. After the fuel sample has burnt for about 20 minutes, the flame is 
extinguished with nitrogen in less than 20 seconds. The burnt sample sheets is cut along the 
longitudinal center line, its surface smoothed, and its thickness measured at the positions along 
the center line where the thermocouples were located. The regression rate is calculated by divid- 
ing the slab thickness difference before and after burning by the effective burning time, which is 
obtained from the temperature history by subtracting the time of pyrolysis front arrival from the 
total burning time. All the experiments were conducted with air as oxidizer, and for flow veloci- 


ties and turbulence intensities ranging from 1 to 4 m/sec and 1 to 20% respectively. 


EXPERIMENTAL RESULTS 


A characteristic example of the measured surface regression rate along the PMMA sheet 
length is presented in Fig. 2 for an air flow velocity of 2 m/s and for several turbulent intensities. 
The regression rate measurements at other flow velocities have similar functional dependences 
as those shown in the figure. From the results of Fig. 2, it is seen that for all turbulence intensi- 
ties the regression rate decreases with the distance from the upstream edge of the PMMA sam- 
ple, the effect being more pronounced as the turbulence intensity is increased. This trend is in 
qualitative agreement with the theoretical predictions from flat plate boundary layer mass burn- 
ing analyses!.7.8, at least for laminar flow. In this process, heat transfer from the flame to the 
solid determines the rate of solid pyrolysis, and consequently of mass burning, and since the 
boundary layer grows thicker downstream, the flame moves away from the surface and the heat 


transferred by convection to the solid is reduced downstream. As a consequence, the surface 


4) 


regression and the mass buming rates decrease as the distance from the PMMA sheet leading 


edge increases, in agreement with the results of Fig. 2. 


The results of Fig. 2 and all the other measurements, also show that the flow turbulence 
affects the regression rate substantially, and that the regression rate can more than double when 
the turbulence intensity is increased from 1% to 20%. Qualitative information about the tur- 
bulence effect can be obtained from Schlieren images of the reacting zone and the thermal boun- 
dary layer. Two Schlieren images are shown in Fig. 3 to illustrate the characteristic differences 
between the laminar flow and the flow with 15% turbulence intensity. In the laminar flow, the 
flame is relatively smooth and stable, and moves away from the surface with the distance from 
the upstream edge. As the flow turbulence is raised to 15%, the flame becomes much more 
unstable and the turbulent vortexes curl the flame and bring it closer to the fuel surface. From 
the photographs, it can be inferred that more heat is transferred to the fuel surface at higher tur- 
bulence intensities because of the proximity of the reaction zone to the fuel surface, and that this 


is the primary factor in the observed regression rate increase with the turbulent intensity. 


The air flow velocity effect on the regression rate along the PMMA sheet is presented in 
Fig. 4 for a flow turbulence intensity of 10%. It is seen from this figure that throughout the fuel 
length, the regression rate is an increasing function of the flow velocity. Similar results are also 
obtained at other flow turbulence intensities. The functional dependence of the regression rate 
on the flow velocity is also in agreement with the predictions from boundary layer mass burning 
analyses!.7.8. As the flow velocity is increased, the boundary layer becomes thinner and the 
flame moves closer to the fuel surface. The result is an enhancement of the convective heat 
transfer from the flame to the fuel and therefore is an increase of the regression rate, as obtained 


in the measurements of Fig. 4. 


DATA CORRELATION 


The qualitative agreement between the regression rate measurements and the predictions of 
the boundary layer analyses suggests the correlation of the experimental data in terms of the 


non-dimensional variables that can be extracted from these analyses. In Fig. 5, the regression 
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rate data is correlated in terms of a non-dimensional burning rate m”/p u.. and the Reynolds 
number, as proposed by Emmons!. The data correlated are for all the flow velocities tested and 
turbulence intensities of 1, 10 and 20%. The results of Fig. 5 show that for a given turbulence 
intensity, the regression rate data correlate well with the selected non-dimensional parameters, 
which suggests that the turbulence intensity should be included as one of the non-dimensional 
parameters in the correlation. In this regard, the analysis of the Smith! for the convective heat 
transfer in a turbulent stagnation point flow, provides a potential method to incorporate the tur- 
bulence intensity in an analysis of the present mass burning problem and to deduce the proper 
non-dimensional parameter to correlate the experimental data. Following Smith’s analysis, an 
eddy kinematic viscosity and diffusivity is defined in the form €=a(u7/u)..u. y, so that it 
includes the turbulence intensity and the shear stress, and heat and mass fluxes are then 
expressed in terms of the laminar, plus the turbulent eddy contribution. With the coordinate on- 
gin at the upstream fuel sheet edge, the flat plate boundary layer conservation equations that 


describe the problem can then be wnitten in the form 


putt + pdt 2/pivee = (2) 
pul «pv = Bl piare z| (3) 


where Z is the mixture fraction? formulated in terms of the Shvab-Zeldovich temperature and 
species coupling functions and where infinite chemistry is assumed. ’:22 The boundary conditions 
for the above equations are that atx =0, y 9 ©, u=u., v=0, and Z=0, and that at y= 0, u=0, 
and Z= 1. The mass burning rate is given by m” =p, r= pv = (q. + q,")/L, where the convec- 


tive surface heat flux is given by q.” =— A (Brp — Bro)(AZ/Ay)y=0- 


Introducing the similarity variable n = yVu../vx, and assuming that there is a locally similar 


solution of the problem such that the stream function can be expressed as y= Vv u. x f(n) and 
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the mixture fraction as Z = Z (n), Eqs. (1) to (3) are reduced to 


[ + too Re,/2 " f+ oe Re, f” + 5 ffa=() (4) 


[ ; : 
e + Gres Rel j Zz’ + E ft+a (as ce Pisa 0) (5) 


with boundary conditions f(0)=f(0)=0, f(e)=1, Z(O)=1, and Z(e%o)=0. The non- 


dimensional mass burning rate 1s then given by 


Peery = Ret? FE ho - 6) 


From Eqs. (4) to (6), it is seen that the solution of the problem has as parameter the product 
(u’/u).. Re,/*, and therefore that the value of (0Z/dn)n=0, and consequently of the mass burning 
rate, will be a function of this flow parameter. Thus, it can be inferred that (u’/u).. Re// should 
be the proper parameter to correlate the experimental data since it describes the flow characteris- 
tics through the turbulence intensity and Reynolds number, and appears as parameter in the solu- 
tion of the governing equations of the problem. 

Correlation of the laminar flow surface regression rate data in terms of the non-dimensional 
mass burning rate (Eq. (6)) and the square root of the Reynolds number is presented in Fig. 6. It 
is seen that there is a linear dependence between the two non-dimensional parameters as 
predicted in laminar, flat plate, boundary layer models of the problem!-’.8. From Fig. 6, the fol- 


lowing expression is deduced for the mass burning rate 
et = 0.28 Re}? (7) 
which is the same expression that is obtained in non-reactive laminar flows for the surface 
Nusselt number?3. 
Correlation of the turbulent flow surface regression rate data according to the predictions of 
Eq. (6), gives an approximately linear dependence between the non-dimensional mass burning 


rate and the square root of the turbulent flow parameter, i.e. ((u’/u).. ReJ->)°>, although the data 
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is somewhat scattered. The scatter of the data is reduced if a 0.8 power is used in the Reynolds 
number, rather than the 0.5 power predicted in the analysis. This is reasonable since in non- 
reacting turbulent flows, it is also found that the Nusselt number correlates with the 0.8 power of 
the Reynolds number?3. The results of the correlation of the turbulent flow surface regression 
rate measurements in terms of the non-dimensional mass burning rate (Eq. (6)) and the new tur- 
bulent flow parameter ((u’/u).. Re2:8)9 > is presented in Fig. 7 for all tests conducted. It is seen 
that all the experimental data correlates very well with these non-dimensional parameters and 
that there is a linear relationship between them. From Fig. 7, the following explicit expression is 


deduced for the mas burning rate 


SIH , 0.5 
mse os 0.8 
Met 4.5 [ y= Rex | (8) 


which provides the relationship between the mass burning rate and the fuel and oxidizing flow 
properties, and corroborates the potential predictive capabilities of the theoretical model intro- 


duced above. 


The correlation of the data in terms of convective heat transfer parameters only indicates 
that for these particular experiments, radiation heat transfer from the flame to the solid has a 
secondary role. It should be noted, however, that the scale of these experiments is small and that 
as the size of the fuel sheet is increased, flame radiation will become increasingly impor- 
tant!2-14. The model introduced above could be extended to approximately include flame radia- 
tion effects by applying the approach followed by Delitchatsios24 in his analysis of the burning 


of large vertical walls in turbulent, natural convection flows. 


CONCLUSION 

Conducting mass burning rate experiments with non-charring, homogeneously burning 
solid fuel in a forced flow with grid induced turbulence, has permitted us to determine the effect 
of the flow parameters on the mass burning rate, and to obtain information useful in the develop- 


ment of theoretical models of solid fuel burning in turbulent flows. The results of the experi- 


45 


ments show that flow turbulence can have a strong effect on the mass burning rate and that this 
effect is due primarily to the enhancement of the convective heat transfer at the fuel surface 
caused by the flame fluctuations. Correlation of the data shows that this effect can be expressed 
in terms of a turbulent flow parameter that includes both the Reynolds number and the tur- 
bulence intensity. It also shows that boundary layer models that consider infinite chemistry and 
incorporate the flow turbulence through an eddy viscosity and diffusivity may be used to predict 
mass burning rates, at least in turbulent flows with a turbulent structure similar to that generated 


with grids. 
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NOMENCLATURE 


A Pre-exponential factor in Arrhenius rate 
a Stagnation point flow stretch, a = du._/dx 
b Parameter in solid heating equations, b = 2h/(rkpc)!/2 
c Specific heat of the solid 
C; Constante1 =12253 
Cp Specific heat of the gas 
D Diffusion coefficient 
d Solid fuel thickness 
Ie Activation energy 
F Fuel 
g Gravity 
h Heat transfer coefficient 
h, Linearized radiation heat transfer coefficient 
I Inert 
k Thermal conductivity 
L Heat of vaporization 
lr Length of the solid heated region ahead of the pyrolysis front 
To Vaporization rate 
n Reaction order 
O Oxidizer 
Products 
Pp Pressure 
Q Heat of combustion 
Oen Heat release rate 
a Heat flux 
as External radiant flux 
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Ub 


Uso 


Heat flux from the flame to the solid 
Cnitical heat flux for pilot ignition 
Radiation from the flame to the solid 
Reradiation from the surface 
Surface heat flux 

Universal gas constant 

Radius 

Temperature 

Flame temperature 

Frozen flow temperature 

Ignition temperature 

Initial solid temperature 

Pyrolysis temperature 

Surface temperature 

Gas initial temperature 

Time 

Ignition time 

Gas induction time 

Solid pyrolysis time 

Velocity parallel to the solid surface 
Free flow velocity 

Forced flow free stream velocity 
Turbulent velocity fluctuations 
Flame spread rate 

Velocity normal to the solid surface 
Molecular weight of specie i 


Reaction rate 


6] 


x Coordinate parallel to the solid surface 

NG Mass fraction of specie i 

y Coordinate normal to the solid surface 

Greek Symbols 

(04 Non-dimensional parameter in gas phase ignition analysis, a = E/RT.. 


Non-dimensional parameter in gas phase ignition analysis, B = T.. cp/YFs Q 


fF Non-dimensional Rae aua es in gas phase ignition analysis, 
T = 4c70(2 — B)/e*(1 — B)? 

) Boundary layer thickness 

€ Small expansion parameter in gas phase ignition analysis, 
E= Rie Ear 

¢ Stretched normal coordinate in gas phase ignition analysis, 
¢ = (1/2e)erf [(2apcp/K,)!* y/(2(1 — exp(—4at))]. 

Jats Characteristic Damkohler number in the stagnation point flow ignition analysis, 
Eqa@aii)): 

xy Characteristic Damkohler number in the flat plate flow ignition analysis, 
Eqs(2.26): 

MT Viscosity 

Vv Kinematic viscosity 

V; Stoichiometric coefficient 

fe) Density 

Oo Non-dimensional parameter in gas phase ignition analysis, 


o = A,(1 — exp(—4at))/(2c2a). 


y Perturbation of the temperature in gas phase ignition analysis, Eq. (2.14) 
Subscripts 

b Buoyant 

cr Critical value 

F Fuel 

f Flame 

fr Frozen flow 
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S 


co 


no-subscript 


Gas phase 

Inert 

Specie 1 

Ignition 

Oxidizer 

Initial solid condition 
Products 

Pyrolysis or vaporization 
Surface 

Initial gas condition 


Solid phase 
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The flaming ignition of a combustible material, and the subsequent spread of flames over its 
surface and steady burning, are critical processes in the initial development of a fire. They deter- 
mine the ease of fire initiation in a given scenario, the materials that will be involved initially in 
the fire and how fast they will become involved, the growth of the fire, and the rate of heat 


release. 


A fundamental understanding of these processes is necessary to develop material flamma- 
bility tests methods that are meaningful and that provide reliable information about the behavior 
of the material in an actual fire. It is also necessary to cevelop the models that are needed to 
predict the fire behavior of the material under different conditions and as input in overall models 
of the development of fires in a variety of scenarios. It is also necessary to produce engineering 
formulas that, although simplified, could be used to predict the fire behavior of materials and the 


development of an actual fire with reasonable accuracy. 


For these reasons, the objective of this chapter is to first provide a coherent description of 
the mechanisms controlling the ignition, flame spread and mass burning of a solid combustible, 
so that a fundamental understanding of these process can be achieved. Afterward, theoretical 
models of these individual processes are developed which include the dominant controlling 
mechanisms, but that are simplified enough to provide explicit formulas that can then be used in 


practical applications and overall models of actual fire development. 


Although wind or air currents induced by openings in the building can drive a fire by forced 
convection, most fires are driven by natural convection flows. Unfortunately it is difficult to test 
in natural convection the effect of the environment on the fire behavior of the materials because 
it requires large scale testing and the control of large gas volumes. For this reason most of the 
detailed studies of the effect of the environment on the ignition and spread of fire over solid 
materials have been conducted under forced flow conditions. Since these studies are very help- 
ful in understanding the mechanisms controlling the development of a fire, and the mechanisms 
are basically the same for both flow conditions, the forced flow studies are presented first so that 


a basic understanding of the specific process under study can be developed. Afterward, their 
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results are extended or modified to describe the more common case of a naturally induced fire. 


2.1 IGNITION OF A SOLID COMBUSTIBLE 


There are different ways that a solid combustible can be ignited and consequently different 
mechanisms that will be influential in the solid ignition process. The ignition, or combustion 
reaction initiation of a combustible material can occur either in the solid or the gas phases. In 
fires, the later case (flaming ignition) is the most important since it may lead to the spread of the 
fire. The gas phase ignition of a solid combustible is generally the combined result of an exter- 
nally imposed heat flux (radiation and convection) that causes the gasification of the solid, and 
the presence of conditions (in the gas or external to it) that will lead to the onset of a sustained 
combustion reaction between the vaporized fuel and the oxidizer gas. If the reaction is initiated 
by an ignition source (open flame, electrical spark, flying ember, etc.) the ignition is normally 
referred to as piloted ignition. If ignition occurs without a pilot the process is normally referred 
to as spontaneous, or auto, ignition. For this type of ignition to occur the gaseous oxidizer/fuel 
mixture must be at elevated temperature. Gas phase ignition caused by a high surface tempera- 
ture or by a surface reaction, such as smoldering or char oxidation, can also be included in this 
Category. 

Several reviews of the ignition of solid combustibles have been published describing the 
different experimental observations and theoretical models of the process (Kanury, 1971,1987, 
Williams, 1985, Drysdale, 1985, Vilyunov and Zarko, 1989), and the reader is referred to these 
works for additional information about the problem. The objective in this work is provide a 
coherent description of the mechanisms controlling the ignition of solid combustibles with an 
emphasis on the development of fires, and to introduce some recent results that may contribute to 
our present knowledge of the subject. In this regard we will first discuss the experimental results 
of Niioka et. al. (1981), which we feel show very well the basic mechanisms that control the 
flaming ignition of a solid combustible. Next we will present a simplified theoretical model of 
the problem which is simple enough to permit deriving explicit formulas describing these 


mechanisms and apply it to other ignition cases. Finally we will discuss the corrections 
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necessary to extrapolate these formulas to more realistic and complicated cases. 


2.1.1 IGNITION CONTROLLING MECHANISMS 


Niioka et. al. (1981) studied experimentally the spontaneous ignition of several solid poly- 
mers in a high temperature oxidizing gas flow. In the experiments, the flat surface of a hemis- 
pherical solid sample was suddenly exposed to a hot, stagnation point, gas flow (Fig. 2.3a), and 
the ignition delay time was measured for different values of the incoming flow velocity (stretch 
rate) and the gas temperature and oxygen concentration. The initial surface heat flux was calcu- 
lated from this data using stagnation point flow boundary layer analysis. The results for the igni- 
tion delay times as a function of the initial surface heat flux, or the stretch rate, are shown in Fig. 
2.1. It is seen that for all fuels tested the ignition time first decreases, reaches a minimum and 
then increases as the flow velocity is increased. For a given fuel, the minimum ignition time 
decreases and occurs at larger flow velocities, as the gas temperature and oxygen concentration 


are increased. 


These results suggest that two primary, competing, mechanisms control the solid fuel igni- 
tion process (Niioka et. al., 1981). One is the heating and gasification of the solid, and the other 
the onset of the gas phase chemical reaction. For example, analyzing the effect of the flow velo- 
city on these mechanisms; as the flow velocity is increased the surface heat flux increases and 
consequently the time required to heat up and pyrolyze the solid decreases. Thus, the influence 
of this mechanism can be graphically represented by the descending line in Fig. 2.2 (pyrolysis 
time). In the gas phase, on the other hand, the increase in the gas velocity results in a decrease 
of the flow residence time (here, the inverse of the stretch rate), which, if comparable with the 
chemical time, will delay the onset of the chemical reaction and eventually prevent its initiation. 
Thus, the effect of this mechanism can be graphically represented by the ascending line in Fig. 
2.2 (induction time). The solid (total) ignition delay time is the sum of the solid pyrolysis time, 
and the gas induction time, as it is sketched in Fig. 2.2. If the gas phase chemical time is very 
short compared with the flow time, the induction time will be also very small and the solid igni- 


tion delay will be controlled by the heat transfer to the solid. This will occur for low flow 
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velocities and/or high temperatures and oxygen concentrations (descending branch in ignition 
lines of Fig. 2.1). For high velocities and/or low temperatures and oxygen concentrations, how- 
ever, the solid ignition delay will be dominated by the onset of the gas phase chemical reaction 


(ascending branch in ignition lines of Fig. 2.1). 


The above concepts can be extended to other geometrical configurations and other sources 
of ignition. In fact, these concepts have been stated before by other investigators although 
perhaps in a less evident or graphical form than that reported in the work of Niioka et. al. (1981). 
A configuration often used in experiments is the flat plate boundary layer type flow (Fig. 2.3), 
where the surface heat flux and the ratio of the flow to the chemical time decrease with the dis- 
tance from the plate leading edge because of the thickening of the boundary layer, thus the igni- 
tion time and the location along the piate of ignition will depend on the flow conditions. For 
conditions when the gas induction time is very short, ignition will be controlled by heat transfer 
to the solid and will occur near the plate’s leading edge. If the induction time is increased 
(increasing the flow velocity or its oxygen concentration for example) the ignition of the solid 
will be further delayed and will occur further downstream from the plate leading edge. The 
presence of a pilot flame will also affect the ignition process by locally reducing the gas induc- 
tion time. The solid heating will consequently be the dominant controlling mechanism of solid 
ignition process, unless the flow velocity is very large or the oxygen concentration very low. 


Some aspects of these different problems are discussed below. 


2.1.2 AUTO IGNITION 


2.1.2.1 High Temperature Stagnation Point Flow 


In this section a simplified model describing the above ignition controlling mechanisms is 
developed. Although the model describes specifically the ignition problem of Niioka et. al. 
(1981), it incorporates all the mechanisms controlling the solid ignition process and can be 
applied to other more complicated ignition problems if appropriated modifications are made. It 
has been selected first for discussion because it is one of the easiest problems to model and there 


is experimental data for comparison with the theory. The model includes a number of 
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simplifying assumptions so that explicit formulas for the solid heating and gas induction times 
can be obtained. The effect of other more complicated, and more realistic, aspects of the prob- 


lem on the analysis predictions is briefly discussed. 


The simplified ignition model to be developed is schematically shown in Fig. 2.3a. A solid 
fuel slab at a constant initial temperature is suddenly exposed in one of its surfaces to a hot, stag- 
nation point flow, oxidizing gas i.e., to an external constant heat flux. The heat flux is absorbed 
at the solid surface, heats the solid and eventually causes it to pyrolyze. The pyrolyzed fuel 
vapor convects and diffuses outwards, mixes with the hot oxidizer gas and forms a combustible 
mixture near the solids surface. The high gas temperature favors the initiation of a gas phase 
combustion reaction between the fuel vapor and oxidizer that, if strong enough to overcome the 
potential heat losses to the surrounding will result in the generation of a sustained combustion 
reaction (flame) over the combustible surface. 

In this simple model two distinct processes are identified; one is the heating and pyrolysis 
of the solid, and the other the onset of the gas phase reaction. With further simplification, the 


two processes can be analyzed independently as follows. 


a) Solid heating and gasification 


Assuming that the solid heating process is uni-dimensional, that the properties are constant, 
that the solid pyrolysis occurs at the surface and is described by a zeroth order Arrhenius law, 
and that there is no surface oxidation or material charring, the governing solid phase equations 


are 


2 
pe aE =k . (2.1) 
rh” = A p exp(—par) (2.2) 


The initial and boundary conditions for Eq. (2.1) are at t=0, T=T, at y=0,k a FRAP = 


Olee ny 


q;” and at y =d, : 
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Since the activation energy in Eq. (2.2) is large for most materials, the pyrolysis rate is 
small for surface temperatures lower than a certain value, Tp, and it can be assumed that the 
solid pyrolysis takes place at a constant surface temperature Tp, which is defined here as the 
"pyrolysis" temperature. The solution of the above equations is further simplified if the surface 


heat flux is assumed constant and if the solid is considered either thermally thick or thin. 


Thermally thick solid: 


In this case the solution of Eq. (1) and its boundary conditions gives the following expres- 
sion for the time required for the surface temperature to reach the pyrolysis temperature 


(Carslaw and Jaeger, 1959), i.e. the pyrolysis time. 


tk pc (Tp- Tp)? 


lb = ———___f,— G23) 


4q,” 
Thermally thin solid: 


A lumped formulation solution of Eq. (1) gives for the pyrolysis time 


ey DAA (2.4) 


This case is less important because in most practical situations, the fuel behaves as thick (Quin- 


tiere, 1984), and for brevity it will be given less attention in this work. 


Equations (2.3) or (2.4) give the descending line in Fig. 2.2 labeled pyrolysis time accord- 
ing to whether the solid is thermally thick or thin. Solids of intermediate thicknesses will have 
curves that lay in between the above two limiting cases. If the gas phase induction time is very 
small (low velocity and high oxygen concentration, and high gas temperature or piloted ignition) 
the pyrolysis time will be nearly equal to the solid ignition time, and Eqs. (2.3) and (2.4) can be 
solely used to predict ignition times as a function of the process parameters. Because the 
fuel/oxidizer mixture becomes flammable almost immediately after solid pyrolysis starts, the 


pyrolysis temperature is often defined as the "ignition" temperature. Although such a definition 
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is not physically correct (Alvares and Martin, 1971), it can be very useful in some practical 


applications, yielding fairly accurate results. 


In the particular case of the experiments of Niioka et. al. (1981) the solid fuel is thick and 
the surface heat flux is uniform across the surface because the gas flow is a stagnation point flow 
but 1s time dependent because the surface temperature varies with time. Approximating the sur- 
face heat flux calculations by assuming a constant surface temperature equal to the pyrolysis 
temperature, then Eq. (2.3) will describe approximately the solid heating process with the sur- 


face heat flux given by (Arpaci and Larsen, 1984) 


Gs” =C1 (Pg cp kg a)? (T..- Tp) - o € (T - TS) (2.5) 


where it has been assumed that the gas does not absorb radiation and that the surrounding walls 


are at T,. Substituting Eq. (2.5) in (2.3), gives for the thick solid pyrolysis time 


t= rs eer (2.6) 

C1 (Pg Cp ky a)! (T..- Tp) - 0 € (Tp _T8) 
which predicts, at least qualitatively, the experimental results of Fig. 2.1 for low stretch rates, 
that is, the initial part of the descending branch of the ignition line. Unfortunately most of the 
data was taken for conditions where both solid and gas phase mechanisms are important, and 
thus the comparison of Eq. (2.6) and the experimental results can be only qualitative. The pyro- 


lysis time for a thin fuel is obtained by substituting Eq. (2.5) into (2.4). 


If both the convective and re-radiation heat fluxes are of the same order of magnitude the 
surface may reach thermal equilibrium before attaining the pyrolysis temperature in which case 
ignition will not occur. This condition determines a minimum, or critical, heat flux for ignition 


given by equating the denominator of Eq. (2.6) to zero. 


b) Gas phase reaction 


If the ambient oxygen concentration is low or the gas velocity is high, the gas induction 


time may be significant and has to be considered in the prediction of the solid ignition times. 
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The mechanisms leading to the gas phase ignition under these circumstances are complicated 
and difficult to predict. Phenomenologically the process is as follows: after solid pyrolysis, the 
fuel vapor (pyrolyzate) leaves the surface and is diffused and convected outwards mixing with 
the ambient oxidizer and creating a flammable mixture near the solid surface. The composition 
of this mixture varies in time and space depending on the gas thermophysical properties, the 
geometry and the characteristics of the flow. If the mixture temperature is increased, either by 
heat transfer from the hot ambient gas, a pilot or any other mechanism, the combustion reaction 
between the fuel vapor and oxidizer gas may become strong enough to overcome the heat losses 


to the solid and ambient and become self-sustained at which point flaming ignition will occur. 


The theoretical modeling of this process includes the solution of the mass, momentum, 
energy and species equations in conjunction with the chemical rate equations, which in general 
is a formidable task. For example, the laminar flow formulation for the stagnation point, convec- 
tive, solid ignition problem studied by Niioka (1981) which is one of the simplest geometries 


possible, is 


Mass 
oe + P22) + a =( (251) 
Momentum 
Energy 
Species 
as ui eve 1 go pps +Vi Ww (2.10) 


7A 


Chemistry 
N N 
> Vix Wi- SVT is m= mi (2.11) 


i=1 i=1 


With boundary conditions 
att=oandy 9, u=u.,V=Vo,T,=To, Yi= Yio,Yr=0 


aC y =O eee ~ (ST) + ri” La ky (GB), T,=Ts, m”=p,v 


Pv ¥i=P,D Si, i#zF, pav(l — Yr)= pp D SE 
The above problem, as formulated, requires a numerical solution which is not easy to 
obtain. A major difficulty in the above formulation is the correct description of the gas phase 
chemistry, and the subsequent solution of the coupled transport and chemical equations. A com- 
mon practice to simplify the solution of this type of problem is the assumption of a single step 
overall reaction with a reaction rate of the Arrhenius type in which case Eq. (2.11) is replaced by 


Vp WF + Vo Wo > Vp Wp (2:12) 
with 


w =A, p" YF Yo exp (-E/RT) (2.13) 


This assumption results in solutions that emphasize the transport aspects of the process and lim- 
its the predictive capabilities of the models. Recently, with the availability of powerful comput- 
ers, multi-step reactions are being incorporated to the models in an attempt to predicts the 


combustion reactions more accurately. The task is still formidable and the solutions limited. 


The introduction of the method of matched asymptotic expansions for large activation ener- 
gies into the solution of combustion problems (Linan, 1974, Williams, 1985) has permitted the 
derivation of explicit expressions for the critical conditions for ignition and the delay time. This 
method makes use of the fact that a large number of fuels (particularly those with a hydrocarbon 


‘ base) have large activation energies, and that because of the Arrhenius character of their 
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combustion reaction the reaction rate is very small for temperatures below a critical value. As a 
consequence the solution of the problem can be divided into two solutions coupled in space and 
time; one non-reactive (frozen or outer solution) that applies to the low temperature regions in 
the mixture, and another reactive (inner solution) that applies to a very small region of the mix- 
ture where the temperatures are highest. Although the resulting solutions are only approximate, 
primarily because of the simple chemistry assumed in the models, they still provide important 


information about the effect of the problem parameters on the ignition process. 


The frozen solution is obtained from the non-reactive form of Eqs. (2.7)-(2.13) and pro- 
vides information about the transport of heat and mass to the reaction zone. Its role on the final 
solution of the problem is primarily the specification of the problem’s Damkohler number and 
the boundary conditions to match the inner and outer solutions. It is also necessary to calculate 
the temperature and species distributions if a second order approximation of the solution is 


needed. 


The reactive solution is the most important one, and provides information about the onset 
of the combustion reaction. Since the reaction rate is very sensitive to the temperature, the reac- 
tion will occur in a small region where the temperatures are highest and the problem can be 
reduced to analyzing the variation of the perturbation of the temperature, y, around the higher 
values of the frozen temperature. The temperature and species concentrations in the inner region 


can then be written in the form (Niioka, 1981) 


T=Ty+eT pW (2.14) 


and 
Yi = Yig — (R T)/Q) ey (2.15) 


where the frozen temperature and species concentrations Ts, and Yj, are in first order equal to 
T.. and Yj... and where € is a small parameter that appears in the expansion of the Arrhenius 
exponent and is of the form € = R T2/T p E. The concept behind this perturbation analysis is that 


if the perturbed temperature increases rapidly the reaction rate will also increase rapidly insuring 


i 


that ignition will occur. Since the changes in the perturbed temperature will occur in a very 
small region, the space coordinate must be expanded by defining a new stretched variable using 


€ as the normalization parameter. The stretched space coordinate is, in this case, 
= (1/26 )-erk | y 2ap,cp/kg y/2((1 — exp (—4a ») . The resulting problem formulation 1s reduced 
to the following energy equation with the perturbed temperature as dependent variable 
Dade 
a Se + (S—B w) exp (y—¢)=0 (2.16) 


with boundary conditions y(0) = d w(ee)/d € = 0. The condition at e is obtained from matching 
the inner and outer solutions. In the above equation, 6 is a non-dimensional time given by 
oO = A, (1 — exp (—4at)/(cz ©) where Ag is the characteristic Damkohler number of the problem 


given by 
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The solution of Eq. (2.16) gives for the induction time (Niioka, 1981) 


=|] In af (2.18) 


where I’ is given by F=4c3 4 [2 — B)/e2(1 — 62) and a= E/RT.. and B=T..c,/Yrs Q. The 
above parameters are further detailed in the nomenclature. Equation (2.18) gives the ascending 
line in Fig. 2.2 labeled gas induction time. For small values of the flow stretch the ignition time 
becomes independent of a and has a minimum given by t,,=I/(4a Aj.) and as the stretch 
increases the ignition time increases and becomes infinite at the value of the stretch at which 


Ag =I. This last condition also determines the critical Damkohler number for ignition 


Noor = T° | (2.19) 


If the Damkohler number is below this critical value, ignition will not occur. From Eqs. (2.17) 
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and (2.19), it is seen that ignition may not occur due to a too low gas temperature (major factor), 
low oxygen concentration, or large stretch. This critical Damkohler number for ignition is often 


derived from steady state analysis of the problem (Fernandez-Pello and Law, 1982). 


For large Damkohler numbers (A, >>I) the gas phase induction time becomes inversely 
proportional to the Damkohler number, t;, = '/(4 a A.) and becomes negligible for large enough 
Ay. This will occur mainly for high gas temperatures and moderate values of the flow velocity 
and oxygen concentration. The presence of a pilot flame could be sufficient for this limit to be 
reached, and in fact this is the case of the pilot ignition experimental studies described below. 
Although the presence of the pilot practically ensures that ignition will occur because the reac- 
tion rate is extremely sensitive to temperature and only very little fuel vapor is needed to initiate 
the reaction (lower flammability limit), it should be noted that if the oxygen concentration is too 
low or the gas velocity too high the presence of a pilot may not be enough to cause the ignition 


of the gas mixture. 


The dependence of the induction time on the gas properties can be deduced from Eq. (2.18) 
and the definition of the Damkohler number (Eq. 2.17). The model predictions can be compared 
at least qualitatively with the experimental results of Breillat et. al. (1978) and Breillat and Van- 
telon (1980a, 1980b) who, using a still gas experimental apparatus, measured the induction time 
of small polymer samples as a function of the ambient gas temperature, oxygen concentration 
and pressure. The experimental results show that the induction time decreases as the gas tem- 
perature, oxygen concentration and pressure are increased in qualitative agreement with the 


predictions of Eq. (2.18). An example of these measurements is shown in Fig. 2.4. 


c) Ignition delay time 
The solid’s ignition delay time is given by the sum of the solid pyrolysis time (Eq. (2.6)) 
and the gas induction time (Eq.(2.18)), i.e., 


nk pc(Tp- To)? 1 wey: 3 
ie= a a In [ x (2.20) 
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This expression describes the "ignition" line in Fig. 2.2, and predicts the existence of a minimum 
of the solid equation delay time with the flow stretch, as was observed experimentally and dep- 
icted in Fig. 2.1. The predicted dependences on the gas temperature and oxygen concentration 
also agree qualitatively with the results of Niioka et. al. (1981). It is also seen that Eq. (2.20) 
predicts two limiting, or critical, heat fluxes (or stretch) for ignition. A lower limit that is given 
by the condition that the solid must gasify, and an upper limit that is given by the condition that 


the Damkohler number must be larger than a certain value. 


The above simplified analysis has been developed with the objective of extracting the basic 
mechanisms that control the solid ignition process, and to provide explicit formulas for the solid 
heating and ignition time. More accurate analysis can be developed by numerically solving the 
coupled solid and gas phase governing equations. An example of such an analysis is the recent 
work of Wang and Yang (1990) for the ignition of a PMMA sphere in a forced oxidizer flow. 
The work also includes experimental measurements of the ignition delay time as a function of 
the gas velocity and oxygen concentration. Both theoretical and experimental results confirm 
the existence of a minimum of the solid ignition time with the flow velocity and the decrease of 
the ignition time as the ambient oxygen concentration is increased, in qualitative agreement with 


the predictions of Eq. (2.20). 


2.1.2.2 High Temperature Flat Plate Flow 


One aspect of the above analysis that must be kept in mind is that the surface heat flux may 
vary with the spatial coordinates. For example in flat plates the convective heat transfer 
coefficient (free or forced) decreases with the distance from the plate leading edge. Thus if con- 
vection is the mechanism of solid heating, ignition will occur first near the plate leading edge. 
However, if convection plays a cooling role as it happens in a plate heated by radiation, the 
solid’s temperature will be higher downstream from the leading edge, and the value of ignition 
time and critical heat flux will depend on the fuel sample size or position of the pilot. The same 
applies if there is a transition from laminar to turbulent flow somewhere along the plate. Thus 


care should be taken when extrapolating small scale data to predict the large scale material 
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ignition behavior. 


If the gas phase transport parameters vary with the location along the solid surface, the 
above analysis and the overall solid ignition problem becomes more complicated because of the 
added spatial effects. For example, in the case of the ignition of a solid by a hot gas in a flat 
plate flow (Fig. 2.3b), the solid pyrolysis will start first at the plate’s leading edge where the sur- 
face heat flux is larger. However, the ignition of the fuel vapor-oxygen mixture may not occur 
there because the gas velocity is largest in this region (lowest residence time), but somewhere 
downstream from the leading edge where as the boundary layer thickens, the velocity gradients 
and the heat losses become smaller, and it is easier for the reaction to take place. In fact ignition 
may not occur because the plate is not long enough to sustain the reactions. Thus, there are addi- 
tional limiting conditions in this case that are determined by the space variables. The variation of 
the ignition location along the flat plate with the flow parameters has been observed experimen- 
tally by Kashiwagi et. al. (1971). They observed that the ignition location moved downstream 
from the plate leading edge as the flow velocity was increased or as the oxygen concentration 
decreased. This behavior has also been predicted theoretically in the steady state, large activa- 


tion energy asymptotics analysis of Fernandez-Pello and Trevino (1985). 


No transient analyses of the problem have been published to date, although an approximate 
solution can be obtained by extending the analysis of Niioka (1981) to the flat plate flow case as 
follows (August and Fernandez-Pello, 1990). Neglecting heat conduction in the x direction and 
approximating the surface heat flux with a constant value by taking the surface temperature 
equal to T,, the solution of the solid heating problem is then given by Eq. (2.3). For laminar 


flow, the surface heat flux can be approximate by 
qs” = (c3 Pg Cp ky u./x)!/2 (T.. — Tp) —O€s Ty - T?) (2:21) 


and the pyrolysis time becomes 
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which describes the descending solid heating line in Fig. 2.2. It should be noted that in this case, 
increasing x has the same result on the predicted pyrolysis time as decreasing u... This is due to 
the respective effects of u.. and x on the thickness of the thermal boundary layer and conse- 


quently on the surface heat flux. 


The gas phase analysis becomes very similar to that developed in section 2.1.2.1(b) but 


with a characteristic Damkohler number given by 
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So asic, RP ay RTS is! 
and the induction time given by 
tin= Ge in [ x 4 (2.24) 


which describes the ascending gas induction time of Fig. 2. Here x/u.. is the flow residence time 
and understandably have opposite effects on the induction time than on the pyrolysis time 
because the flow residence time must balance the chemical time for the combustion reaction to 
occur. Physically, as the thermal boundary layer grows (x increasing or u.. decreasing) the con- 
vective heat losses from the incipient combustion reaction to the surface and ambient decrease 


helping the onset of the reaction. 


The solid ignition time is obtained by the addition of the solid pyrolysis and gas induction 


times, i.e. 


tig = _ URNS yt - @Inl- ~ (2.25) 

(C3 Pg Cp ky un/x)!? (T.. _ Tp) — 6 €s (T4 - TA) 
It is seen from Eq. 2.25 that in this case, there also are two limiting conditions for ignition. One 
for small Damkohler numbers where the ignition time is controlled by gas phase chemical kinet- 
ics and becomes infinite for A, = T° and another for large Damkohler numbers where the ignition 
time is controlled by the solid heating and pyrolysis and becomes infinite when the surface re- 


radiation balances the convective heat flux. The trends predicted by Eq. (2.25) for the variation 
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of the ignition location with the flow parameters agree qualitatively with the experimental obser- 


vations of Kashiwagi et. al. (1971). 


The above analysis can be readily extended to describe the ignition of the solid in a natural 
convection flat plate flow by simply replacing u.. with a buoyant velocity u, (vg)! or by 


expressing the surface heat flux and Damkohler number in terms of the Grashoff number. 


2.1.2.3 High Surface Temperature 


When the solid is heated by an external source of thermal radiation, the solid surface tem- 
perature may reach a high enough value to cause the gas phase ignition. This type of ignition 
generally occurs with charming materials because the char left after the fuel is pyrolyzed can 
reach very high temperatures due to radiation absorption or char oxidation. If the solid is of the 
charring type, a layer of carbonaceous material is formed at the solid’s surface and can grow 
inward as the more volatle components of the material are pyrolyzed by the external radiant 
flux. This char layer deters the solid pyrolysis by both absorbing the external heat flux and 
deterring the free flow of volatile to the surface, thus modifying the heat of vaporization and the 
surface temperature. The char can also react with the ambient oxidizer in the form of an exoth- 
ermic surface reaction that may result in an increase of the surface temperature. By either radia- 
tion absorption or oxidation, the char can reach a very high temperature and act as an ignition 
source itself. This form of ignition has a particular interest in initiation of wood fires, and nor- 
mally occurs when there is a sudden enhancement of the in-depth solid gasification due to an 
increase in the external heat flux or the oxygen concentration at the surface (with air currents for 
example). 

The ignition model developed above can also be applied to describe this type of ignition, 
although the formation of char will affect the solid heating process. Ignition of the 
vapor/oxidizer gas mixture will occur very near the fuel surface where the temperatures are 
higher. For short ignition times (high surface heat flux, low charring materials) the above model 
should apply if ignition occurs before the material charring becomes important, and Eq. (2.20) 


should describe well the ignition time by simply substituting T;.. by T, in the induction term of 
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the equations. The effect of charring, and the transition from a char surface reaction to flaming, 
however, are complicated processes and it is difficult to incorporate them in the model. A possi- 
ble approach to incorporate charming into the model is to use some of the analyses developed to 


describe mass burning of charring materials (section 2.3). 


A number of studies have been conducted on this type of ignition, most of them with cellu- 
lose sheets exposed to an external radiant flux in natural convection. Kanury (1987) has studied 
and reviewed this problem in great detail and the reader is referred to that work for more details 
on the subject. Particularly interesting are the results of Martin (1965) for the spontaneous igni- 
tion of vertical cellulose sheets under an external radiant flux The investigator plotted the igni- 
tion delay time results in terms of two nondimensional parameters; one that is the ratio of the 
total incident radiation to the energy required to heat up the solid to the pyrolysis temperature, 
and the other that is the ratio of the incident radiant flux to the characteristic solid conductive 
heat flux, and identified the ignition regions controlled by convection cooling, solid conduction 
heating, and solid pyrolysis. The plot, therefore, represents the processes involved in the heating 
and pyrolysis of the solid, and should be described by a version of Eqs. (2.1) and (2.2) that 
would include material charring and natural convection heat transfer. Such an analysis has been 


developed by Kanury (1987) with good qualitative agreement with the expenments. 


The effect of the gaseous environment on the ignition of the cellulose sheets has been stu- 
died experimentally by Alvares and Martin (1971). These experiments complement the work of 
Martin (1965), and show that the gas properties (pressure and oxygen concentration) influence 
the solid ignition temperature and ignition time, thus verifying the role of the gas phase reaction 
on the overall ignition time. The measurements show that the surface temperature at the time of 
ignition, and the ignition time, decrease as the ambient oxygen concentration or the pressure are 
increased, the latter in qualitative agreement with the predictions of Eq. (2.20). Alvares et. al. 
(1970) also studied the influence of the sheet size on the cellulose ignition, and observed that 
once the boundary layer had developed the ignition occurred at the top of the sheet, and the igni- 


tion time decreased as the height of the sheet was increased, also in qualitative agreement with 
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the predictions of Eq. (2.20). 


It should be noted that the model developed above describes this type of ignition only for 
some particular cases, and that the analysis must be extended to account for other effects that 
may appear in the process, such as in-depth absorption of radiation (Ohlemiller and 
Summerfield, 1971), charring (Kanury 1987), surface oxidation (Kashiwagi and Ohlemiller, 
1982 and Kashiwagi et. al. 1985) and thermal stability and melting (Kashiwagi and Omon, 
1988), among others. The analyses of Kashiwagi (1974), Kindelan and Williams (1975,1977), 
and Gandhi and Kanury (1986) are examples of analyses that could be used to extend the present 


model to describe some of the above indicated effects. 


2.1.2.4 Gas Phase Radiation Absorption 


Another aspect of the gas phase ignition process that also deserves attention is the potential 
absorption of radiation by the gas phase. This can occur because the fuel vapor is composed by 
heavy, radiation absorbing, molecules, and/or the presence of soot or other articulates (carried by 
the flow). The absorption of radiation by the gas can affect both the solid heating process, by 
reducing the radiant flux reaching the solid surface, and the gas phase ignition process by 
increasing the temperature of the fuel/oxidizer mixture. This type of ignition has been studied 
experimentally by Mutoh et. al. (1978) and Kashiwagi (1979a, 1979b, 1981) by radiating 
PMMA samples with a CO) laser. The experiments show that the solid eerie rate decreases 
due to the absorption of radiation by the pyrolyzate and that spontaneous ignition of the gas can 
occur if the radiation intensity is above a certain minimum (of the order of 60 W/cm’). 
Although these radiant fluxes are too large to be expected in a normal fire, the importance of the 
observation is that they show that radiation can be absorbed in the gas phase and that it can 
affect the solid ignition process. 

Numerical analyses of the gas phase ignition of a vaporizing solid surface by absorption of 
radiation by the gas have been developed by Amos et. al. (1988), Amos and Fernandez-Pello 
(1988) and Park (1989). The results agree qualitatively with the experimental observations and 


provide information about the minimum radiant intensity for ignition as a function of the 
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problem parameters. Particularly interesting is the prediction of the events leading to the initia- 
tion of ignition and the subsequent establishment of a diffusion flame over the solid surface 
reported by Amos et. al. (1988). The predicted evolution of the species concentration profiles is 
shown in Fig. 2.5. Ignition is predicted to occur first at a certain distance from the surface that 
depends primarily on the absorption coefficient and the radiant flux. Subsequently two premixed 
flames propagate away from the ignition zone, one toward the surface through an increasingly 
rich mixture until it cannot propagate because the mixture is too rich or is quenched by the wall, 
and the other away from the surface through an increasingly leaner mixture unt] the mixture is 
stoichiometric, at which point the flame stabilizes itself in the form of a diffusion flame. This 
series of events, although predicted for this particular problem, should also describe other igni- 


tion problems, particularly those with a pilot as ignition source. 


2.1.3 PILOT IGNITION 
2.1.3.1 External Radiant Flux 


This ignition problem has been given considerable attention by several investigators during 
the last few years because it simulates a form of ignition that occurs often in fires. Here we will 
discuss primarily the work on the subject by Quintiere and co-workers (1981, 1983, 1984), 
which we find instructive in describing the relationship between pilot ignition and flame spread, 
and in identifying the material properties that have a significant role on these processes. The 
work has also resulted in the development of an important method for determining the flamma- 


bility characteristics of combustible materials. 


The specific problem studied by these investigators is the pilot flame ignition of a vertical, 
thick fuel solid exposed to a constant external radiant flux, and the subsequent lateral flame 
spread over the vertical surface that is exposed to a decreasing radiation with the sample length. 
Since ignition is achieved with a pilot flame, the gas induction time is small and the ignition time 
should be very close to the pyrolysis time. Thus, the solution of the problem is given by Eq. 
(2.3), where the pyrolysis time can be referred to as the ignition time and the pyrolysis tempera- 


ture as the. ignition temperature. Linearizing the surface re-radiation and incorporating the 
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resulting radiation heat transfer coefficient with the convective coefficient, and taken the average 


as constant, the following expression is obtained for the pilot ignition time 


__ ukpe(Ti—T,) 
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where the solid initial temperature has been taken equal to the ambient temperature and the sur- 
face temperature equal to the pyrolysis temperature, qe” is the external radiant flux and h incor- 
porates the convective and linearized radiation heat transfer coefficients. If the external radiant 
flux is small, the solid surface may come into thermal equilibrium before the pyrolysis tempera- 
ture is reached in which case ignition will not occur. This condition leads to the determination 


of the minimum heat flux for solid ignition as 


Gps =h (Tig = bg) (2.27) 


Substitution of Eq. (2.27) in (2.26) gives for the ignition time 


Fi 
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Quintiere et. al. (1983,1984), have measured the piloted ignition time for many building materi- 

als as a function of the external radiant flux, and have successfully correlated most of their meas- 


urements with an expression similar to Eq. (2.28). An example of these correlations is given in 


Fig. 2.6 that has been extracted from Quintiere (1981) and Quintiere and Harkleroad (1984). 


Given the success in the correlation of their ignition results, the above investigators pro- 
posed the use of their experimental procedure to determine material property data that would 
characterize the ignition of the material. Also, instead of using Eq. (2.28), they proposed to 
correlate the data with a simpler version of the equation of the form q’o,ig/de” = F(t) with 
F(t) = b Vt for t<t,, and F(t) = 1 for t > tm, which is basically Eq. (2.28) for qe” > de o,ig and b = 
2h/(x k p c)'/2. The correlation of the experimental data is then used to derive b, t,,, and Giote. 


which are used to comparatively classify the ignition characteristics of the materials and as input 
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data in mathematical fire models. 


In the course of their work these investigators also showed that the size and position of the 
pilot flame can affect the measured value of the ignition time and critical heat flux for ignition. 
Unless the pilot flame is small and removed from the solid surface, it can contribute to the heat- 
ing of the solid and as a consequence the measured values of the above parameters will vary 
with the pilot size and location. On the other hand if the pilot is not positioned within the region 
close to the solid surface where the flammable mixture is developed, ignition will not occur. 
Thus care should be taken to assure that the ignition measurements are not pilot dependent. A 
common solution to this problem is to place a small pilot flame just above of the vertical solid 


surface. 


Several investigators have also experimentally determined the inverse relationship between 
the ignition time and the externally applied radiant heat flux, and deduced the critical value of 
the radiant flux for piloted ignition for a variety of materials ranging from cellulose to polymers 
(Simms, 1963, Alvares, 1975, Tewarson and Pion, 1978). These last studies have been reviewed 
in detailed by Drysdale (1985) and the reader is referred to that work for further details. 
Recently Hasegawa et. al. (1986) and Tewarson and Kahn (1988) have conducted similar experi- 
ments in their studies of the ignition of electrical cable insulations. Although with electrical 
cables the conducting core can play an important role as a heat sink, particularly for long expo- 
sure times to radiation or for thin insulating jackets, the ignition characteristics of the cables are 


generally similar the ones described above. 


A critical parameter in the application of the above ignition model is the pyrolysis or igni- 
tion temperature, Tj,. Its value can be directly measured with thin thermocouples embedded at 
the solid’s surface, although the measurement may not be very accurate because it is difficult to 
place the thermocouple exactly at the surface, and it may also absorb or emit radiation dif- 
ferently than the surface itself. Another method is to experimentally determine the critical heat 
flux for ignition with plots like that of Fig. 2.6 and then use Eq. (2.27) to obtain Tj,, which if 


assumed constant can be used in the above equations to calculate the ignition time. This tem- 
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perature, however, depends on several of the problem parameters, such as surface heat flux, oxy- 
gen concentration, charring, etc., and care should be taken when extrapolating values obtained 
under specific conditions to other situations. The fact that the pyrolysis temperature is not con- 
stant can be seen simply through the solution of Eqs. (2.1) and (2.2), and has been shown for dif- 
ferent problems by Fernandez-Pello and Williams (1977a), Altenkirch (1982), and Di Blasi et. 
al. (1989b). Also relevant to this aspect of the problem are the experimental observations of 
Alvares and Martin (1971) and of Kashiwagi et. al. (1988, 1989) that show that the pyrolysis 
temperature varies with the oxygen concentration and ambient pressure. The latter authors also 
observed that the temperature is different when the ambient gas is pure nitrogen or has oxygen, 


indicating the presence of a surface oxidation reaction that influences the solid pyrolysis process. 


2.2 FLAME SPREAD OVER A SOLID COMBUSTIBLE SURFACE 


The spread of a flame over the surface of a solid combustible is a subject of interest in fire 
safety because it influences the initial fire development and rate of heat release. A schematic of 
the spread of a flame over a solid fuel surface is shown in Fig. 2.7 for two characteristic flame 
configurations, opposed (2.7a) and concurrent (2.7b) flame spread. For the flame to spread, 
enough heat must be transferred from the flame to the unburnt material ahead of the flame to 
pyrolyze the solid. The vaporized fuel is then diffused and convected away from the surface, 
mixing with the oxidizer and generating a flammable mixture ahead of the flame leading edge, 
which is then ignited by the flame. The rate of flame spread is therefore determined by the abil- 
ity of the flame to transfer the necessary heat to pyrolyze the solid and to ignite the combustible 


mixture ahead of it. 


The heat transfer from the flame to the unbumt combustible ahead is strongly dependent on 
the shape of the flame, which in turn is dependent on the characteristics of the gas flow. When 
the gas flow, either naturally induced or forced, opposes the direction of spread (Fig. 2.7a), the 
flow keeps the flame close to the surface downstream of the pyrolysis front, deterring the heat 
transfer ahead of the flame. This type of spread is commonly known as opposed flow flame 


spread and is generally slow. In natural convection, it occurs in downward or horizontal flame 
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flames over the PMMA rods and sheets in a still, oxygen enriched, atmosphere show that the 
flame spread rate in micro-gravity is in between the normal gravity downward and upward flame 
spread rates. This result indicates that in downward spread, the buoyant entrainment of air by 
the flame has a deterring effect on its spread probably due to the convective cooling which dom- 


inates over the increased convective flow of oxidizer into the reaction zone. 


Olson et. al. (1988, 1990) have conducted a series of experiments in a drop tower of the 
spread of flames over thin paper sheets in a low velocity opposed gas flow. The results show 
that there is a minimum flow velocity at which extinction is observed and that as the velocity is 
increased, the spread rate first increases, reaches a maximum of approximately 10 cm/sec and 
then decreases. Measurements of the oxygen concentration at which extinction is observed also 
show that there is a minimum ( 14%) at a flow velocity of 10 cm/sec. These results are interest- 
ing because they suggest that a micro-gravity environment may be more dangerous than a nor- 
mal gravity one from the point of view of material flammability. These experiments should be 
described well by Eq. (2.31) since at micro-gravity, the flow is truly a forced one. At low flow 
velocities, however, Eq. (2.31) predicts a spread rate that is independent of the flow velocity. A 
possible explanation is that in these experiments, the paper sheets behave as a thick fuel rather 
than as a thin one because the solid thermal penetration layer is thinner at these low flow veloci- 
ties. If the paper sheets behave as thick, then Eq. (2.30) applies and predicts better the experi- 
ments. The extinction limit observed experimentally at low flow velocities would be caused by 
surface heat losses according to Eq. (2.30) or (2.31). At low flow velocities, the flame stand-off 
distance is larger and the surface heat flux may become equal to the heat losses (Olson et. al. 
1988, Chen, 1989, Bhattacharjee et. al. 1990). Another possible explanation for the observed 
increase of the spread rate with the flow velocity is that at these low flow velocities, oxygen 
depletion near the reaction zone due to diffusive effects may have an important role in the spread 
of the flame (Olson, 1990). One interesting aspect of this is that, if verified, would add a third 


controlling mechanism of flame spread that would be due to reactant diffusion. 
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2.2.2.3 Concurrent Flow 


In the concurrent mode of flame spread, the flame covers the solid element during the heat- 
ing and pyrolysis processes (Fig. 2.7b), and when the pyrolyzed fuel leaves the surface, is 
rapidly ignited by the flame. Thus, in this case, the ignition process can be viewed as a strong 
piloted ignition or one with an ambient gas at very high temperatures, and that is consequently 
controlled by the solid’s heating and pyrolysis. Thus, the flame spread rate will be given by Eq. 


(2.30) or (2.31) with the gas induction term neglected. 


The length of the solid heated region, /;, is better defined in this case since it is closely 
related to the length of the flame (from the pyrolysis front to the flame tip). The hot post- 
combustion gases downstream from the flame tip also contribute to the solid heating, although it 
appears that this contribution is secondary in comparison to that of the flame itself (Loh and 
Fernandez-Pello, 1984). For non-charring materials, the length of the flame has been found to be 
proportional to the length of the solid burning region (pyrolysis length) (Orloff et. al. 1975, Saito 
et. al. 1986, Tu and Quintiere, 1988) in the form /f = cs J} with 0.8 <n < 1.1. Choosing n = 1 for 


simplicity, the laminar flow concurrent flame spread rate is then given by 


rah 4 C5 ((c3 k p U/L ue (TAs Tp) + qr + der) I 4 


for thick solids, and 
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for thin solids. 

The above equations should predict the concurrent flame spread rate well for conditions 
away from extinction. Experimentally it has been observed that extinction generally occurs first 
at the flame upstream leading edge and then propagates downstream until the whole flame is 
extinguished. The limiting conditions for extinction could be approximately determined with 
aetT (2.52). 


Concerning the surface heat flux, Orloff et. al. (1975) in their study of diffusion flames over 
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large vertical PMMA surfaces, have found that convection is important for flames smaller than 
approximately 20 cm and that radiation becomes dominant for larger flames. Thus, as the flame 
spreads, the solid heating will change from convection initially to radiation later on. From Eq. 
(2.36) it is seen that during the initial convective period, the spread rate will be independent of 
[p, linearly proportional to u.. and proportional to Y 2. (through Ty), in agreement with the exper- 
iments presented in Fig. 2.12. The independence of Vs on /, is due to the increase in the flame 
stand off distance that counteract the increase in the flame length (Mao and Fernandez-Pello, 
1981). As the flame spreads and radiation becomes dominant, the spread rate will accelerate 
exponentially according to Eq. (2.36). The experiments of Apte et. al. (1990) for flames spread- 
ing over large horizontal PMMA sheets confirm the acceleration of the flame as it spreads, 


although in the experiments buoyancy had a strong effect on the spread of the flame. 


For thin fuel sheets, there is fuel burn-out at the upstream fuel leading edge and the pyro- 
lysis length is determined by the rate of propagation of the pyrolysis and burn-out fronts. The 
flame will still be proportional to the pyrolysis length, and since this length will not reach a large 
value, convection heat transfer will probably be the major heat transfer mechanism in this case. 
Thus, Eq. (2.37) will contain the convective term only (no external radiation) with J, = xp — Xp, 
where x, and x, are the distances from the initial sheet leading edge to the pyrolysis and bum- 
out fronts respectively. Equation (2.37) predicts a spread rate that initially increases with time as 
I, increases but that eventually could become constant if the velocities of the pyrolysis and 
burn-out fronts become equal. It also predicts an initial dependence of the spread rate on the 
flow velocity which should become weaker as the burn-out rate increases due to the enhance- 
ment of the burning rate. These, and the predicted linear dependence of the spread rate with the 
oxygen concentration agree qualitatively with the experiments of Loh and Fernandez-Pello 
(1985) presented partially in Fig. 2.13, and with the predictions of the numerical analysis of Di 
Blasi et. al. (1989c). 


The effect of charring. Charting of the solid combustible will have a strong influence on 


the concurrent mode of flame spread because of its direct impact on the pyrolysis length. Unfor- 
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tunately, the problem is too complex to develop a simple model of the process that can be easily 
incorporated in the above analysis. Carrier et. al. (1983) and Di Blasi et. al. (1989c) have incor- 
porated in their numerical analyses of the concurrent flame spread the effect of charring but in a 


too simplified form to describe well the actual charring of the material. 


The effect of charring appears primarily through the fuel pyrolysis rate that, as the char 
layer builds up, starts to decrease and may even not occur if the char retains enough heat to 
prevent the pyrolysis of the material underneath. Thus, charring introduces through the pyro- 
lysis rate, and consequently through the pyrolysis length, a time dependent function into the 
problem that is difficult to solve analytically. Saito et. al. (1986) and Quintiere et. al. (1986) 
have developed analyses for upward flame spread that incorporate approximately a time depen- 
dent burning rate and that seems to predict qualitatively the experimental observations. Such an 
approach could also be used in the forced flow analysis although the resulting flame spread rate 
equation would be somewhat complicated. Also, since the result of charring is the a spread of a 
"charming" or burn-out front, it may be possible to describe the spread rate over charring materi- 
als with a thin solid expression of the form of Eq. (2.37) with /, derived from experimental infor- 
mation. The dependences on the oxygen concentration and flow velocities may also be depen- 
dent on the generation of char since the char layer will influence the thermal behavior of the 


solid. 


The effect of flow turbulence. In the concurrent mode of flame spread the effect of flow 
turbulence is particularly important because the spread of the flame is strongly dependent on the 


flame length and surface heat flux, and they are in turn strongly affected by the flow turbulence. 


Zhou and Fernandez-Pello (1990) have recently conducted an experimental study of the 
effect of flow turbulence on the concurrent mode of flame spread that has produced very infor- 
mative and somewhat unexpected results. The measurements of the rate of flame spread over the 
surface of thick PMMA sheets as a function of the velocity and turbulence intensity of a con- 
current forced air flow are presented in Fig. 2.17. They show that flow turbulence has a strong 


influence on the flame spread rate, with the spread rate decreasing as the turbulence intensity 
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increases, the effect being more pronounced for larger flow velocities. This effect is due pn- 
marily to a strong influence of the turbulence intensity on the flame length, that as seen in Fig. 
2.18 decreases as the turbulent intensity is increased, and that counteracts the increase of the sur- 
face heat flux with the turbulent intensity. The above investigators did not determine the cause 
for the decrease of the flame length when the turbulent intensity was increased, but speculated 
that it was quenching of the diffusion flame due to the strong entrainment of cold ambient air 


into the reaction zone and the interaction with the cold wall. 


The flame spread rate should still be predicted in this case by Eq. (2.36) with the surface 
heat flux and /; reflecting the flow turbulence effect. In order to verify the predictive capabilities 
of Eq. (2.36) the above authors substituted the flame spread rate, flame length, and surface heat 
flux data in Eq. (2.36) and plotted the results as a function of the turbulence intensity and 
obtained a good correlation of the measurements, which again verifies the predictive capabilities 
of Eq. (2.36) and corroborates that the shortening of the flame length with the turbulence inten- 
sity is the primary reason for the observed reduction of the flame spread rate as the turbulence 
intensity is increased. These results are important not only because they introduce new and 
unexpected aspects about the flame spread process not previously known, but because they may 
have a significant influence in the application of flame spreads formulas in models of room fire 


development. 


2.2.3 NATURAL CONVECTION FLAME SPREAD 


Flame spread in natural convection is particularly important in the fire safety field because 
most fires are driven by free or mixed convection flows. For this reason, a large number of the 
flame spread investigations made to date have been conducted in natural convection. Unfor- 
tunately, most of these works have been concerned primarily with solid phase effects and have 
not addressed the gas phase contribution, probably because in free convection it is difficult to 
control the ambient conditions. For natural convection flame spread, the concepts and analysis 
developed for forced flow can be readily applied by simply replacing the forced flow parameters 


with natural convection ones, or more simply, replacing the flow velocity u. with a buoyant 
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velocity. Also, here it is convenient to differentiate between opposed and concurrent flows. 


2.2.3.1 Opposed Flow (Downward) 

Downward flame spread over a vertical surface corresponds to this type of flame spread. 
Eqs. (2.30) and (2.31) can be used to predict the downward flame spread rate by substituting the 
forced flow parameters (q,” and 8) with natural convection ones, or simply by replacing the flow 
velocity u.. with a buoyant velocity of the form (de Ris, 1969) 


ie 1B 
- gee nee 


Pg Cp I 


and 6 (x/g)!4. The resulting equation also predicts a flame spread controlled by a balance 
between solid heating and gas phase reaction. Since in most fires the gravity vector and pressure 
remain constant, the latter controlling mechanisms becomes important under vitiated conditions. 
The predictions of Eq. (2.30) and (2.31) with (2.38) agree qualitatively with the experimental 
measurements of Altenkirch et. al. (1980, 1983) for the spread of flames over thin paper and 
thick PMMA sheets at varied gravity, pressure and oxygen mass fractions. They also contain the 
non-dimensional parameters, (non-dimensional flame spread, and Damkohler number) used by 


these above authors to correlate their data. 


The extension of the analysis to horizontal flame spread is more difficult because radiation 
from the flame to the solid surface ahead of the flame is the dominant mechanism of heat transfer 
(Fernandez-Pello et al., 1980) and consequently the solid surface heat flux and length of the 
heated region depend on a number of parameters such as flame shape, emission, soot concentra- 
tion, etc. that are difficult to predict theoretically. Equation (2.29) is still applicable if a proper 
surface heat flux and solid heated length are used. 

External Radiation. In fires, the combustible material is often preheated by radiation from 
an adjacent burning material. Thus, it is of interest to predict the fire behavior of the material 
under the influence of an external radiant flux. The effect of an externally applied radiation on 


the flame spread rate appears primarily through the initial temperature of the solid and can be 
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incorporated in the above analysis through the value of T,. This has been verified by 
Fernandez-Pello (1977b) that correlated the downward spread rate over PMMA sheets exposed 
to different radiant fluxes in terms of (T, — T,)* with T, being the measured surface temperature 


at the instant of flame arrival. 


Quintiere et. al. (1981, 1983, 1984) using the radiant panel test apparatus, have made exten- 
sive measurements of the rate of lateral flame spread over a large number of building materials 
exposed to an external radiant flux. The lateral flame spread is also described by Eq. (2.29) 
although the length /; and heat transfer coefficients are specific to this particular geometry. The 
above authors resolved the problem of determining these unknown parameters by grouping them 
into a major one that is empirically obtained by correlating the flame spread measurements with 
the analysis predictions. Furthermore, they proposed to use this parameter together with those 
derived for solid ignition to classify the flammability characteristics of the materials. The com- 
bined ignition and flame spread analysis of the material is very interesting and since it provides 
information about key material properties regarding fire initiation and spread, it may prove to be 


an important and useful method for testing material flammability. 


Quintiere et. al. (1981, 1983, 1984) proposed to correlate their lateral flame spread rate 


measurements with an equation for the form 


Vi= ipo tye (2.39) 


which is basically Eq. (2.29) with the surface heat flux and heated length incorporated in 6 and 
with T;, = T, and with the gas reaction term neglected. Eq. (2.39) should therefore predict well 
the flame spread rates for conditions where solid heating is the controlling mechanism but will 
start to fail for those conditions where gas phase chemical kinetics become important, such as in 
vitiated ambients. In order to express the flame spread rate in terms of the external radiant flux 
and to unify the ignition and flame spread analyses, Eq. (2.39) is further modified by substituting 
Ti, with qo,ig using Eq. (2.27), and T, with the external radiation q-”(x) using the solution of the 


solid heating problem. The resulting expression for the flame spread rate is 


22 


1 
Vr= Spy PAL 2.40 
ida : ar . 
: C* (q oig ~ Ge (x)F(t)) 


where C = (4/n)!/*/b. The flame spread rate measurements are conducted for long periods of 
time in which case F(t) = 1, and the data is plotted in terms of q,”(x) in the same graph as the 
one used to plot the ignition time data (Fig. 2.6) and Eq. (2.40) is used to derive the value of C, 


which is considered a material property. 


Plotting the ignition and flame spread data in a single graph as that of Fig. 2.6 results in a 
very useful way to analyze the flammability characteristics of the materials. The above investi- 
gators named these graphs "flammability diagrams" and proposed to use them in conjunction 
with Eqs. (2.28) and (2.40) to obtain the ignition and flame spread properties of the materials 
Sani Tig, b and C. These properties can then be used to classify the materials according to 
their ease of ignition and flame spread. These diagrams also help to graphically explain the rela- 
tionship between ignition and flame spread. In the first place, it is seen that both the ignition 
delay time and the flame spread rate tend asymptotically to infinity from opposite sides when the 
external radiant flux approaches the critical heat flux for ignition, which defines a boundary 
between ignition and flame spread. If the external radiant flux is larger than q’,;, enough 
energy reaches the solid to pyrolyze it and if a pilot is present ignition will occur. If the external 
radiant flux is smaller than q’’oj,, ignition will not occur unless enough heat is transferred to the 
surface from another source (a spreading flame) to bring the total surface heat flux to a value 
larger than q”,;,. If the sum of the external heat flux and the heat flux from the flame is larger 
than q”’o,ig, ignition will occur and the ignition delay time will be given by the ignition branch of 
Fig. 2.6 (or Eq. (2.28)) with q.” equal to the sum of the actual external heat flux plus the heat 
flux from the flame. The resulting flame spread rate will be given by the length of the solid sur- 
face heated by the flame divided by the ignition time, which when plotted in Fig. 2.6 in terms of 


the external radiant flux will give a point in the flame spread curve of the figure (Eq. 2.40). 


In an actual fire, a normal sequence of events leading to the ignition and spread of flames 


over a combustible material is that the solid is heated locally by an external heat flux of magni- 
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tude greater than q’’o,ig and the fuel vapor is ignited by a hot spot. The flame, once estab.ished, 
spreads by itself or assisted by an external heat flux. For some matemials, the heat flux from the 
flame to the solid is not enough to pyrolyze the solid and the flame cannot spread unless is 
assisted by an external heat flux. The minimum external heat flux for flame spread is, therefore, 
another important parameter in the flame spread characteristics of the material and is included 
together with the above indicated ignition and flame spread parameters in the flammability 
classification of the material. The reader is referred to Quintiere et. al. (1984) for further discus- 


sion of the application of the flammability diagrams. 


2.2.4.2 Concurrent Flow (Upward) 


Upward flame spread is perhaps the most important mode of flame spread in fire safety 
because it is often present during the development of a fire and is rapid and hazardous. Consid- 
erable work has been conducted in this area during the last decade, both experimental and 
theoretical, and there is an active research effort currently underway. The reader is referred to 
the reviews of Fernandez-Pello and Hirano (1983) and Fernandez-Pello (1984) for earlier works 


on the subject. 


This type of flame spread is basically the same as the concurrent forced flow flame spread 
presented before except that, since the flow is naturally induced by the flame, the flow parame- 
ters are natural convection ones. Thus, the concepts regarding the pilot ignition characteristics 
of the process and the length of the heated region also apply here, and the flame spread rate is 
given by the solid heating terms of Eqs. (2.30) or (2.31) but with a natural convection surface 
heat flux and the length /; equal to the flame height (from the pyrolysis front to the tip). For lam- 
inar flow this gives the following expression for the rate of upward flame spread rate over a thick 


solid fuel surface 


_ 4¢5((C6 ky Pg Cp 8 (Tr — Too)/(xT..))/* (Te — Tp) + Gr” + Ge”)? 
Vr= eae Te ee (2.41) 


which is similar to that given by Sibulkin and Kim (1977). In this case also, natural convection 
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heat transfer will be dominant initially, but as the flame progresses and the size of the flame 
increases, the flame will become turbulent and radiation will start to take over. Both effects will 
affect the heat transfer and flame length characteristics of the process. Equation (2.41) shows 
that in this case the spread rate is acceleratory even under laminar flow conditions, which is in 
agreement with the experimental observations of Fernandez-Pello (1977b) and Annamalai and 
Sibulkin (1979). This is because the increase of the flame stand off distance as the flame spreads 
and the boundary layer thickens does not counteract the corresponding increase of the flame 
length (Mao and Fernandez-Pello, 1981). As the flame becomes turbulent and the heat transfer 
radiation dominated, Eq. (2.41) predicts that the flame spread becomes more accelerative, also in 
qualitative agreement with the measurements of Saito et. al. (1986) for the upward spread of 


flames over thick PMMA sheets. 


For thin materials Eq. (2.31) becomes 


Mecs(caks Peco eS (lr (x be) ipa) Oe eigpg it 
‘a a a (2.42) 


with /, =X, —Xp. Thus the characteristics of the flame spread will also depend here on the rela- 
tive rates of spread of the pyrolysis and burn- out fronts. The experiments of Markstein and de 
Ris (1973) for the upward spread of flames over thin textile sheets show a spread process that is 
similar to the forced flow measurements presented in Fig. 2.13, with a spread rate that is 
acceleratory initially and later becomes constant as the burn-out and pyrolysis fronts reach the 
same spreading velocity. 

In most practical situations the upward spread of fire is turbulent, radiation dominated and 
the material chars, and although the above equations still apply, they must be adapted to account 
for these effects. The random character of these processes makes the problem extremely 
difficult to solve theoretically, and the adaptation of the above equations to predict the rates of 
flame spread must be based on simplifying assumptions or empirical correlations. An active 
research effort has been underway during the last few years and progress is being made to under- 


stand and predict better this type of flame spread. Saito et. al. (1986, 1989) studies of the upward 
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turbulent spread of flames over wood provides interesting information about the effect of char- 
ring on the upward spread of the flame, and show that for charring materials it may be necessary 
to expose the material to an external radiant flux for the flames to spread. If the radiant flux is 
above a certain minimum value the flames spread at increased rates as the radiant flux is 
increased. Their measurements of the upward flame spread rate over thick particle-board sheets 
as a function of the external radiant flux are presented in Fig. 2.19. An interesting aspect of this 
plot is that it is equivalent to the flame spread part of the flammability diagram of Fig. 2.6 and 
corroborates the suggestion of Hasegawa et al. (1986), that similar flammability diagrams could 


be developed for upward flame spread. 


Saito et. al. (1986) also developed a model of flame spread that includes similar concepts to 
those introduced in this work, and that incorporates the effect of charring through the fuel mass 
burning rate, which in turn determines the flame length through the rate of heat release (Deli- 
chatsios, 1984). They also correlated the flame spread rate data in terms of the ratio of the meas- 
ured flame length and the pyrolysis time, thus verifying the validity of Eq. (2.31). Other recent 
studies and models that incorporate the same concepts are those of Quintiere et. al. (1986) and 


Hasemi (1986) 


Delichatsios (1984, 1986, 1988) has also studied this problem and addressed the need for 
data on pyrolysis rates, flame heights, radiation and surface heat fluxes in order to predict flame 
spread rates. He also developed a flame spread model that incorporates the concepts presented 
here, and that includes simplified expressions for the surface heat flux and flame heights and that 
provides an expression for the spread rate that can be used to correlate flame spread data. Deli- 
chatsios (1988) proposes five parameters, extracted from the spread rate expression, that can be 
determined in a small-scale test and that could be applied to characterize the fire hazard charac- 
teristics of the material. Together with the flammability diagrams, these parameters should per- 


mit a good determination of the flammability characteristics of combustible materials. 
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2.3 MASS BURNING AND HEAT RELEASE RATES 


Once the solid is ignited and the flame spreads over its surface, a diffusion flame is esta- 
blished over the solid and sustains its pyrolysis by heat transfer to its surface. The solid pyro- 
lysis, Or mass burning, rate and the subsequent rate of heat release by the combustion of the 
pyrolyzed fuel and the ambient oxidizer, are important parameters in the development of a fire 
and information about their properties is needed as input in room fire models. They determine 
such aspects of the fire development as flame heights (impingement on walls, flame spread, radi- 
ation to adjacent surfaces, etc.), and the room ambient temperature (ceiling layer depth and tem- 
perature, wall temperature distribution etc.). A considerable number of studies have been con- 
ducted on the subject, and the reader is referred to Drysdale (1985), Pagni (1981) and Sibulkin 
(1988) for additional information on these works. 

For a non-charring material the pyrolysis, or mass burning, rate is given by the solution of 
Eqs. (2.1) and (2.2) and their boundary conditions. If the pyrolysis temperature is taken as con- 
Stant (section 2.1), and it is assumed that the heat flux to the solid interior is small in comparison 
to the heat needed to pyrolyze the solid, then the solution of the problem is considerably 


simplified and the pyrolysis rate will be given by 


m” = (qr a Ges tae Gn’ VL (2.43) 


where the surface heat flux has been separated into the heat transferred from the flame to the 
solid by convection and radiation, an external source of radiation, and the surface re-radiation. 
Assuming that all the pyrolyzed fuel burns with the ambient oxidizer, the total instantaneous rate 
of heat release per unit width will be 
Xp(t) 
Qin= { m”’Qdx (2.44) 
Xp(t) 


where x,(t) and x,(t) will be given by their respective spread rates. 


The above statements concerning flame turbulence and radiation, and material charring also 


97 


apply here, and in fact they also make this problem difficult to solve. An active research effort 
has been underway for the last few years in an attempt to solve some of the different aspects of 
this process and to provide practical information on the rate of burning and heat release for use 
in room fire models. Delichatsios (1984, 1985, 1988) has studied turbulent wall fires and 
developed simplified expressions for flame heights, and mass burning and heat release rates that 
are easily applicable to predict upward flame spread rates and total heat release rates. The 
expressions are similar to those of Eqs. (2.44) and (2.45), and include approximate relations for 
the different surface heat fluxes. Mitler (1988) also developed a model of surface pyrolysis for 
non-charring materials that incorporates heat convection, flame radiation and surface reradiation 
and that can readily be used in room fire models. Kulkarni and Kim (1988) have studied the 
effect on the pyrolysis rate of the heat conducted to the solid interior for thick PMMA slabs and 
concluded that this heat flux is not negligible and that should be considered when calculating 
mass burning rates. Recently Zhou and Fernandez-Pello (1990) have studied the effect of grid 
induced turbulence on the mass burning rates of PMMA sheets in a flat plate air flow, and found 
that the mass burning rate measurements can be correlated in terms of nondimensional parame- 
ters deduced from turbulent boundary layer theory, as it is shown in Fig. 2.20. From the correla- 


tion of the results, the following expressions are deduced for the mass burning rate 


a) laminar flow 


keGh eee) 
xX 


mn” = 0.28 8)’ Re} (2.45) 


and 


b) turbulent flow 


= , 0.5 
i= 4st ete ess ci | (2.46) 


which are in agreement with Eq. (2.43) with the surface heat flux obtained from boundary layer 


diffusion flame analysis (Pagni, 1981). This result is intersting because verifies the applicability 
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of boundary layer analysis to predict flame spread and mass burning rates. The above forced 
flow results can be extended to free convection flow by replacing the flow velocity u. with a 
buoyant flow velocity (Eq. 2.38) or using a surface heat flux based on natural convection boun- 


dary layer analysis. 


The effect of material charring is also very important here since it affects directly the pyro- 
lysis rate, generally making it a time decaying process. A considerable amount of experimental 
work has been conducted on the subject and the reader is referred to Drysdale (1985) and the 
recent review of Atreya (1984) for details on these works. Models have also been developed of 
the pyrolysis of charring materials, both numerically (Kung 1972, Parker 1985) and analytically 
(Kanury 1972, Miller and Ramohalli 1986, and Wichman and Atreya 1988), that address some 
of the aspects of the problem. The process however is very complex, and although some of the 
analyses provide explicit formulas for the pyrolysis rate and char growth that are directly appli- 
cable, there is still a need for further studies on the problem. Another aspect of the problem that 
deserves attention is the potential thermal degradation and oxidation of the material by the 
ambient oxygen, which may affect the value of the pyrolysis activation energy obtained by ther- 
mal degradation measurements in inert atmospheres (Hirata et. al. 1985). Finally, it should be 
kept in mind that in fires, the gas flows are often mixed convective (free and forced) and that the 
formulation of the above different processes may need to be adapted to include such types of 
flows. The analysis of Mao et. al. (1984) is an example of how such a formulation can be imple- 


mented. 


2.4 SUMMARY 


The controlling mechanisms of ignition, flame spread and mass burning, are well under- 
stood and it is possible to derive rate expressions for these processes that can predict well rela- 
tively simple cases (laminar flow, non-charring, homogeneous materials, non-radiating flames, 
etc.). More difficult is to predict these processes under more complex, and more realistic, situa- 
tions were turbulence, material properties, and flame radiation play a role, and there is a need for 


more advances in these areas. For practical applications, however, the simplified expressions 
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can be used to predict complex situations by either approximating the affected parameters or by 


using correlations derived from laboratory tests. 
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